For R beginners

New chunk Ctrl+Alt+I

Execute chunk Ctrl+Shift+Enter

Execute all chunks Ctrl+Alt+R

HTML preview Ctrl+Shift+K

Library preparations

library(readr)
library(dplyr)

Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union
library(tidyverse)
── Attaching core tidyverse packages ────────────────────────────────────────────────────────────────────────────────────────────────────── tidyverse 2.0.0 ──
✔ forcats   1.0.0     ✔ stringr   1.5.1
✔ ggplot2   3.5.1     ✔ tibble    3.2.1
✔ lubridate 1.9.3     ✔ tidyr     1.3.1
✔ purrr     1.0.2     ── Conflicts ──────────────────────────────────────────────────────────────────────────────────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag()    masks stats::lag()
ℹ Use the ]8;;http://conflicted.r-lib.org/conflicted package]8;; to force all conflicts to become errors
library(ggplot2)
library(reshape2)

Attaching package: ‘reshape2’

The following object is masked from ‘package:tidyr’:

    smiths
library(stats)
library(tm)
Loading required package: NLP

Attaching package: ‘NLP’

The following object is masked from ‘package:ggplot2’:

    annotate
library(text2vec)
Registered S3 method overwritten by 'data.table':
  method           from
  print.data.table     
library(textstem)
Loading required package: koRpus.lang.en
Loading required package: koRpus
Loading required package: sylly
For information on available language packages for 'koRpus', run

  available.koRpus.lang()

and see ?install.koRpus.lang()


Attaching package: ‘koRpus’

The following object is masked from ‘package:tm’:

    readTagged

The following object is masked from ‘package:readr’:

    tokenize
library(syuzhet)

Data Import

data_posts <- read.csv("~/4year/2semester/dtII/CSVs/HEIs.csv",
                 colClasses = c(tweet_id = "character"))

# Modifying created_at type so that attribute can be used more easily 
data_posts$created_at <- as.POSIXct(data_posts$created_at,
                              format= "%Y-%m-%dT%H:%M:%S", tz="UTC")

#View(data)
summary(data_posts)
      id              tweet_id             text               type           bookmark_count    favorite_count     retweet_count      reply_count      
 Length:11728       Length:11728       Length:11728       Length:11728       Min.   :  0.000   Min.   :    0.00   Min.   :   0.00   Min.   :   0.000  
 Class :character   Class :character   Class :character   Class :character   1st Qu.:  0.000   1st Qu.:    7.00   1st Qu.:   2.00   1st Qu.:   0.000  
 Mode  :character   Mode  :character   Mode  :character   Mode  :character   Median :  0.000   Median :   20.00   Median :   5.00   Median :   1.000  
                                                                             Mean   :  1.543   Mean   :   60.67   Mean   :  10.62   Mean   :   3.888  
                                                                             3rd Qu.:  1.000   3rd Qu.:   57.00   3rd Qu.:  11.00   3rd Qu.:   3.000  
                                                                             Max.   :418.000   Max.   :41655.00   Max.   :4214.00   Max.   :2317.000  
                                                                                                                                                      
   view_count        created_at                       hashtags             urls            media_type         media_urls       
 Min.   :      5   Min.   :2022-08-01 03:05:11.00   Length:11728       Length:11728       Length:11728       Length:11728      
 1st Qu.:   2643   1st Qu.:2022-10-19 12:56:27.00   Class :character   Class :character   Class :character   Class :character  
 Median :   6240   Median :2023-01-29 08:26:30.00   Mode  :character   Mode  :character   Mode  :character   Mode  :character  
 Mean   :  14182   Mean   :2023-01-30 07:39:34.96                                                                              
 3rd Qu.:  16058   3rd Qu.:2023-05-05 14:16:43.25                                                                              
 Max.   :7604544   Max.   :2023-08-31 20:50:01.00                                                                              
 NA's   :4840                                                                                                                  

Data Pre-Processing

# Count of how many entries each HEI has
number_posts <- data_posts %>%
              group_by(id) %>% summarise(count = n())

number_posts

Since complutense only has 1 entry we can’t learn anything from it, so we removed it

data_posts <- data_posts[data_posts$id != "complutense.csv", ]

Visualization of all posts, just tweets and just replies

number_posts <- data_posts %>%
              group_by(id) %>% summarise(posts = n())

number_tweets <- data_posts[data_posts$type == "Tweet", ] %>%
              group_by(id) %>% summarise(tweets = n())

number_replies <- data_posts[data_posts$type == "Reply", ] %>%
              group_by(id) %>% summarise(replies = n())

print(number_posts)
print(number_tweets)
print(number_replies)

Calculating the percentage of tweets and replies based on all posts

# Merging the counts of tweets and replies with the count of posts
data_ratio <- merge(number_posts, number_tweets, by = "id", all = TRUE)
data_ratio <- merge(data_ratio, number_replies, by = "id", all = TRUE)


data_ratio$percentage_tweets <- round(((data_ratio$tweets / data_ratio$posts) * 100), 2)
data_ratio$percentage_replies <- round(((data_ratio$replies / data_ratio$posts) * 100), 2)

data_ratio <- data_ratio[, c("id", "percentage_tweets", "percentage_replies")]

data_ratio$percentage_replies[is.na(data_ratio$percentage_replies)] <- 0

print(data_ratio)

NA removal

Function to visualize the number of NAs in all columns

na_count <- function(){
  # Counting the number of NA values for each column
  na_count <- colSums(is.na(data_posts))
  
  # Creating a new data frame with the NA counts
  na_counts_table <- data.frame(Column = names(na_count), NA_Count = na_count)
  
  print(na_counts_table)
}

Calculations of view, favourite, retweet and reply percentiles and visualization of NAs in all columns

data_posts <- data_posts %>%
  group_by(id) %>%
  mutate(view_percentile = ntile(view_count, 100),
         favorite_percentile = ntile(favorite_count, 100),
         retweet_percentile = ntile(retweet_count, 100),
         reply_percentile = ntile(reply_count, 100)) %>%
  rowwise() %>%
  mutate(avg_percentile = round(mean(c(view_percentile, favorite_percentile, retweet_percentile, reply_percentile), na.rm = TRUE), 2))

na_count()

data_percentile <- data_posts[, c("id", "view_percentile", "favorite_percentile", "retweet_percentile", "reply_percentile", "avg_percentile")]

print(data_percentile)

Calculation of the maximum number of views for each HEI

max_view_counts <- tapply(data_posts$view_count, data_posts$id, max, na.rm = TRUE)

print(max_view_counts)
      duke.csv       epfl.csv        goe.csv    harvard.csv  leicester.csv manchester.csv        mit.csv         sb.csv   stanford.csv    trinity.csv 
        307969         105095          15455        2982704          47838         317086        7604544         607498         222593         205333 
        wv.csv       yale.csv 
        109265         143108 

Removal of NAs

# From view count
data_posts$view_count <- ifelse(
  is.na(data_posts$view_count),
  round(max_view_counts[data_posts$id] * (data_posts$avg_percentile / 100)),
  data_posts$view_count)

# From view percentile
data_posts$view_percentile <- ifelse(
  is.na(data_posts$view_percentile),
  data_posts$avg_percentile,
  data_posts$view_percentile)

Visualization of NAs in all columns

na_count()

Data Processing

Function to calculate average posts

average_posts <- function(timeframe){
  # Calculation of the timeframe between earliest and latest post for each HEI
  date_range <- data_posts %>%
    group_by(id) %>%
    summarise(min_date = min(created_at),
              max_date = max(created_at)) %>%
    mutate(num_days = as.numeric(difftime(max_date, min_date, units = timeframe)))
  
  # Naming the column respecting the timeframe
  column_name <- paste0("avg_posts_per_", timeframe)
  
  # Calculation of the number of posts per day for each HEI
  posts_per_timeframe <- number_posts %>%
    left_join(date_range, by = "id") %>%
    mutate(!!column_name := round((posts / num_days), 2))
  
  print(posts_per_timeframe)
  return(posts_per_timeframe)
}
posts_per_day <- average_posts("days")
posts_per_week <- average_posts("weeks")

Plot for the average number of posts per day for each HEI

barplot(posts_per_day$avg_posts_per_days,
        names.arg = posts_per_day$id,
        main = "Average Posts per Day",
        xlab = "HEI",
        ylab = "Average Number of Posts",
        ylim = c(0, max(posts_per_day$avg_posts_per_days) + 1),
        las = 2,
        col = "#3498DB")

# Adding text labels over each bar and aligning it with the center of each bar 
text(x = barplot(posts_per_day$avg_posts_per_days, plot = FALSE),
     y = posts_per_day$avg_posts_per_days,
     labels = round(posts_per_day$avg_posts_per_days, 2),
     pos = 3)

Plot for the average number of posts per week for each HEI

barplot(posts_per_week$avg_posts_per_weeks,
        names.arg = posts_per_week$id,
        main = "Average Posts per Week",
        xlab = "HEI",
        ylab = "Average Number of Posts",
        ylim = c(0, max(posts_per_week$avg_posts_per_weeks) + 5),
        las = 2,
        col = "#E74C3C")

text(x = barplot(posts_per_week$avg_posts_per_weeks, plot = FALSE),
     y = posts_per_week$avg_posts_per_weeks,
     labels = round(posts_per_week$avg_posts_per_weeks, 2),
     pos = 3)

Defining the intervals of time for the academic year

intervals <- list(
  interval1 = as.POSIXct(c("2022-08-31", "2022-12-15")),
  interval2 = as.POSIXct(c("2023-01-04", "2023-04-01")),
  interval3 = as.POSIXct(c("2023-04-14", "2023-06-15"))
)

Function to check if a date falls within a given interval of time and apply appropriate Boolean

check_interval <- function(date) {
  for (i in 1:length(intervals)) {
    interval_start <- intervals[[i]][1]
    interval_end <- intervals[[i]][2]
    if (date >= interval_start & date <= interval_end) {
      return(TRUE)
    }
  }
  return(FALSE)
}
data_posts$academic_year <- sapply(data_posts$created_at, check_interval)
print(data.frame(id = data_posts$id, academic_year = data_posts$academic_year))

Function to count number of posts and average per day during academic time and vacation time

analyze_posts <- function(academic_year_filter) {
  # Filtering the data based on the academic_year_filter
  filtered_data <- data_posts %>%
    filter(academic_year == academic_year_filter)
  
  # Count of days for each HEI
  unique_days <- filtered_data %>%
    group_by(id) %>%
    summarise(unique_days = n_distinct(as.Date(created_at)))
  
  # Count of posts for each HEI
  number_posts_boolean <- filtered_data %>%
    group_by(id) %>%
    summarise(count = n())
  
  # Naming the column respecting the time period
  time <- ifelse(academic_year_filter, "academic_time", "vacation_time")
  column_name <- paste0("avg_posts_in_", time)
  
  # Combination of data and calculation of average posts per day
  combined_data <- left_join(unique_days, number_posts_boolean, by = "id")
  combined_data <- combined_data %>%
    mutate(!!column_name := round((count / unique_days), 2))
  
  print(combined_data)
  return(combined_data)
}
data_posts_academic <- analyze_posts(TRUE)
data_posts_vacations <- analyze_posts(FALSE)

Plot for the average number of posts during academic time for each HEI

barplot(data_posts_academic$avg_posts_in_academic_time,
        names.arg = data_posts_academic$id,
        main = "Average Posts during Academic Time",
        xlab = "HEI",
        ylab = "Average Number of Posts",
        ylim = c(0, max(data_posts_academic$avg_posts_in_academic_time) + 5),
        las = 2,
        col = "#34495E")

text(x = barplot(data_posts_academic$avg_posts_in_academic_time, plot = FALSE),
     y = data_posts_academic$avg_posts_in_academic_time,
     labels = round(data_posts_academic$avg_posts_in_academic_time, 2),
     pos = 3)

Plot for the average number of posts during vacation time for each HEI

barplot(data_posts_vacations$avg_posts_in_vacation_time,
        names.arg = data_posts_vacations$id,
        main = "Average Posts during Vacation Time",
        xlab = "HEI",
        ylab = "Average Number of Posts",
        ylim = c(0, max(data_posts_vacations$avg_posts_in_vacation_time) + 5),
        las = 2,
        col = "#D35400")

text(x = barplot(data_posts_vacations$avg_posts_in_vacation_time, plot = FALSE),
     y = data_posts_vacations$avg_posts_in_vacation_time,
     labels = round(data_posts_vacations$avg_posts_in_vacation_time, 2),
     pos = 3)

Data preparation for dates

# Creating new table that contains a new column for the day of the week
data_posts_days <- data_posts %>%
  mutate(day_of_week = weekdays(created_at))

# Selecting only the id, created_at, and day_of_week columns for the new table
data_posts_days <- data_posts_days %>%
  select(id, created_at, day_of_week)

# Create column hour from created_at
data_posts_days$created_hour <- as.numeric(format(data_posts_days$created_at, "%H"))

print(data_posts_days)
# Grouping by id and day_of_week, then counting the number of posts
number_posts_days <- data_posts_days %>%
  group_by(id, day_of_week) %>%
  summarise(count = n())
`summarise()` has grouped output by 'id'. You can override using the `.groups` argument.
# Grouping by id, day_of_week and day created at, then counting th enumber of tweets
number_posts_per_day <- data_posts_days %>%
    mutate(created_date = as.Date(created_at)) %>% 
    group_by(id, day_of_week, created_date) %>%
    summarize(count = n())
`summarise()` has grouped output by 'id', 'day_of_week'. You can override using the `.groups` argument.
# Finding for each HEI the average count of posts per day
average_number_posts_per_day <- number_posts_per_day %>%
  group_by(id, day_of_week) %>%
  summarise(average_count = round(mean(count), 2))
`summarise()` has grouped output by 'id'. You can override using the `.groups` argument.
print(number_posts_days)

Highest and lowest posts

# Finding the HEI with the lowest count of posts per day
lowest_count <- number_posts_days %>%
  group_by(day_of_week) %>%
  slice_min(order_by = count) %>%
  select(day_of_week, id, count)

# Finding the HEI with the highest count of posts per day
highest_count <- number_posts_days %>%
  group_by(day_of_week) %>%
  slice_max(order_by = count) %>%
  select(day_of_week, id, count)

# Combine the results
high_low_HEI <- bind_rows(lowest_count, highest_count) %>%
  arrange(day_of_week)

print(high_low_HEI)

Plot for the highest and lowest count of posts per day for each day of the week

ggplot(high_low_HEI, aes(x = day_of_week, y = count, fill = id)) +
  geom_bar(stat = "identity", position = "dodge") +
  geom_text(aes(label = count),
            position = position_dodge(width = 0.9),
            vjust = -0.5,
            size = 3) +
  labs(title = "Highest and Lowest Count of Posts per Day for Each Day of the Week",
       x = "Day of the Week", y = "Count") +
  scale_fill_manual(values = rainbow(length(unique(high_low_HEI$id)))) +
  theme_minimal() +
  theme(legend.title = element_blank())

Average of posts

# Finding the HEI with lowest and highest averaged count of posts per day
high_low_average_HEIs <- average_number_posts_per_day %>%
  group_by(day_of_week) %>%
  filter(average_count == max(average_count) | average_count == min(average_count)) %>%
  arrange(day_of_week, ifelse(average_count == min(average_count), average_count, -average_count))

print(high_low_average_HEIs)

Plot for the highest and lowest average count of posts per day for each day of the week

ggplot(high_low_average_HEIs, aes(x = day_of_week, y = average_count, fill = id)) +
  geom_bar(stat = "identity", position = "dodge") +
  geom_text(aes(label = round(average_count, 2)),
            position = position_dodge(width = 0.7),
            vjust = -0.5,
            size = 3) +
  labs(title = "Highest and Lowest Average Count of Posts per Day for Each Day of the Week",
       x = "Day of the Week", y = "Average Count") +
  scale_fill_manual(values = rainbow(length(unique(high_low_HEI$id)))) +
  theme_minimal() +
  theme(legend.title = element_blank())

Favourite hour and day

favourite_day_hei <- number_posts_days %>%
  group_by(id) %>%
  top_n(1, count) %>%
  arrange(id)

print(favourite_day_hei)
# Calculating the number of posts per hour
number_posts_hours <- data_posts_days %>%
  group_by(id, created_hour) %>%
  summarise(count = n()) %>%
  ungroup()
`summarise()` has grouped output by 'id'. You can override using the `.groups` argument.
# Identifying the favorite hour for each HEI
favourite_hour_hei <- number_posts_hours %>%
  group_by(id) %>%
  top_n(1, count) %>%
  arrange(id)

# Adding new columns to categorize the hours and assign numerical value for clustering
favourite_hour_hei <- favourite_hour_hei %>%
  mutate(time_of_day = case_when(
    created_hour >= 5 & created_hour < 12 ~ "Morning",
    created_hour >= 12 & created_hour < 18 ~ "Afternoon",
    TRUE ~ "Night"
  ),
  time_of_day_value = case_when(
    created_hour >= 5 & created_hour < 12 ~ 1,
    created_hour >= 12 & created_hour < 18 ~ 3,
    TRUE ~ 5
  ))

print(favourite_hour_hei)

Heatmaps

Function to plot heatmap for various HEIs

heatmap_maker <- function(target_id){
  # Filtering data for the specific HEI
  target_data <- data_posts_days %>%
    filter(id == target_id)
  
  # Grouping by day of the week and hour, and counting the number of tweets
  tweet_counts <- target_data %>%
    group_by(day_of_week, created_hour) %>%
    summarise(num_posts = n())
  
  # Plotting heatmap
  ggplot(tweet_counts, aes(x = day_of_week, y = created_hour, fill = num_posts)) +
    geom_tile() +
    scale_fill_gradient(low = "white", high = "blue") +
    labs(title = paste("Post Heatmap for", target_id),
         x = "Day of the week",
         y = "Hour of the day")
}

Plot of heatmap for each HEI

heatmap_maker("duke.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("epfl.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("goe.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("harvard.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("leicester.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("manchester.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("mit.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("sb.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("stanford.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("trinity.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("wv.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

heatmap_maker("yale.csv")
`summarise()` has grouped output by 'day_of_week'. You can override using the `.groups` argument.

Hashtags

# Transforming empty strings into NA
data_posts$hashtags[data_posts$hashtags == ""] <- NA

# Table with number of unique hashtags and percentage of usage
hashtags <- data_posts %>%
                group_by(id) %>%
                summarise(count = n(),
                          na = sum(is.na(hashtags)),
                          unique_hashtags = length(unique(hashtags)),
                          hashtag_percentage = round(((count - na) / count * 100), 2))

print(hashtags)

Plot for the count of unique hashtags for each HEI

barplot(hashtags$unique_hashtags,
        names.arg = hashtags$id,
        main = "Unique Hashtags for Each HEI",
        xlab = "HEI",
        ylab = "Count of Unique Hashtags",
        ylim = c(0, max(hashtags$unique_hashtags) + 50),
        las = 2,
        col= "#16A085")

text(x = barplot(hashtags$unique_hashtags, plot = FALSE),
     y = hashtags$unique_hashtags,
     labels = round(hashtags$unique_hashtags, 2),
     pos = 3)

Plot for the usage of hashtag for each HEI

barplot(hashtags$hashtag_percentage,
        names.arg = hashtags$id,
        main = "Hashtags Percentage for Each HEI",
        xlab = "HEI",
        ylab = "Hashtags Percentage",
        ylim = c(0, max(hashtags$hashtag_percentage) + 30),
        las = 2,
        col= "#F1C40F")

text(x = barplot(hashtags$hashtag_percentage, plot = FALSE),
     y = hashtags$hashtag_percentage,
     labels = round(hashtags$hashtag_percentage, 2),
     pos = 3)

URL usage

# Transforming empty strings into NA
data_posts$urls[data_posts$urls == ""] <- NA

# Table with number of post, number of NA and url percentage of usage
url_usage <- data_posts %>%
                group_by(id) %>%
                summarise(count = n(),
                          na = sum(is.na(urls)),
                          url_percentage = round(((count - na) / count * 100), 2))

print(url_usage)

Plot for the usage of hashtag for each HEI

barplot(url_usage$url_percentage,
        names.arg = url_usage$id,
        main = "Urls Percentage for Each HEI",
        xlab = "HEI",
        ylab = "Urls Percentage",
        ylim = c(0, max(url_usage$url_percentage) + 10),
        las = 2,
        col= "#8E44AD")

text(x = barplot(url_usage$url_percentage, plot = FALSE),
     y = url_usage$url_percentage,
     labels = round(url_usage$url_percentage, 2),
     pos = 3)

Text

data_posts_content <- data_posts %>%
            select(id, text)

# Counting number of words
data_posts_content <- data_posts_content %>%
  mutate(num_words = lengths(strsplit(text, "\\s+")))

# Grouping by HEI and calculate average, minimum, and maximum values of number of words
data_posts_content_metrics <- data_posts_content %>%
  group_by(id) %>%
  summarise(average_num_words = mean(num_words),
            min_num_words = min(num_words),
            max_num_words = max(num_words))
print(data_posts_content_metrics)

Plot for the average, maximum and minimum values of words for each HEI

ggplot(data_posts_content_metrics, aes(x = id, y = average_num_words)) +
  geom_point(aes(color = "Average")) +
  geom_errorbar(aes(ymin = min_num_words, ymax = max_num_words, color = "Range"), width = 0.2) +
  scale_color_manual(values = c("Average" = "#1976D2", "Range" = "#EF5350")) +
  labs(title = "Word Count Summary by HEI",
       x = "HEI",
       y = "Number of Words",
       color = "Metric") +
  theme_minimal() +
  theme(axis.text.x = element_text(angle = 90, vjust = 0.5, hjust = 1))

Posts Classification

Function to cleanup text

cleanup <- function(docs, spec.words=NULL){
  # lowercase
  docs <- tm_map(docs, content_transformer(tolower))
  # rm numbers
  docs <- tm_map(docs, removeNumbers)
  # rm english common stopWords
  docs <- tm_map(docs, removeWords, stopwords("english"))
  # if stopwords are specified as a character vector
  if(!is.null(spec.words))
    docs <- tm_map(docs, removeWords, spec.words)
  # rm punctuations
  docs <- tm_map(docs, removePunctuation)
  # rm extra white spaces
  docs <- tm_map(docs, stripWhitespace)
  # lemmatizing text
  docs <- tm_map(docs, lemmatize_words)
  
  docs
}

Fucntion to pair words

pairing <- function(word_dict, paired_words) {
  if (nrow(word_dict) != length(paired_words)) {
    stop("The number of rows in word_dict and the length of paired_words should be equal.")
  }
  
  word_dict <- tibble::tibble(word = word_dict$value)
  result <- tibble::tibble(word = word_dict$word, category = paired_words)
  return(result)
}

Function to classify text

classify_text <- function(text, word_pair) {
  # Tokenizing the text and converting to lowercase
  words <- tolower(unlist(strsplit(text, "\\W+")))
  
  # Finding the frequent words in the text
  freq_words <- words[words %in% word_pair$word]
  
  if (length(freq_words) == 0) {
    return("Unknown")  # Returning Unknown if no frequent words are found
  }
  
  # Getting the corresponding categories
  categories <- word_pair$category[word_pair$word %in% freq_words]
  
  # Sorting categories by frequency in descending order and return the most frequent one
  return(names(sort(table(categories), decreasing = TRUE))[1])
}

Function to create corpus

corpus_maker <- function(text) {
  texts <- text$text
  vc <- VectorSource(texts)
  corpus <- Corpus(vc)
  
  return(corpus)
}

Function to find frequent terms in texts

freq_terms <- function(clean_text, number) {
  dtm <- DocumentTermMatrix(clean_text)
  dtm.tfidf <- weightTfIdf(dtm)
  
  mdf <- as_tibble(as.matrix(dtm.tfidf))
  
  mdf.freq <- mdf %>%
  select(findFreqTerms(dtm, number)) %>%
  summarise_all(sum) %>%
  gather() %>%
  arrange(desc(value))

  mdf.freq$key <- 
    factor(mdf.freq$key,
           levels = mdf.freq$key[order(mdf.freq$value)])
  
  word_dictionary <- as_tibble(mdf.freq$key)
  
  ggplot(mdf.freq, aes(x=key, y=value)) +
    geom_bar(stat="identity") +
    labs(x="terms", y="freq") + coord_flip()
  
  return(word_dictionary)
}

Duke

duke_text <- subset(data_posts_content, id == "duke.csv")
duke_corpus <- corpus_maker(duke_text)

duke_clean_text <-cleanup(duke_corpus, c("new", "will", "change", "north", "can", "first", "year", "carolina", "years", "summer", "meet", "one"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
duke_word_dictionary <- freq_terms(duke_clean_text, 30)

duke_word_category <- c("Image", "Education", "Education", "Image", "Image", "Research", "Research", "Image", "Image", "Engagement", "Research", "Research", "Image", "Education", "Education", "Society", "Education", "Engagement")

duke_word_pair <- pairing(duke_word_dictionary, duke_word_category)

print(duke_word_pair)

# Applying the classify_text function to each text
duke_text <- duke_text %>%
  mutate(category = sapply(text, classify_text, word_pair = duke_word_pair))

print(duke_text)

epfl

epfl_text <- subset(data_posts_content, id == "epfl.csv")
epfl_corpus <- corpus_maker(epfl_text)

epfl_clean_text <- cleanup(epfl_corpus, c("new", "amp", "can", "will", "one", "now", "–", "via", "portrait"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
epfl_word_dictionary <- freq_terms(epfl_clean_text, 30)

epfl_word_category <- c("Image", "Research", "Image", "Research", "Research", "Research", "Education", "Image", "Research", "Research", "Education", "Research", NA, "Education", "Education")

epfl_word_pair <- pairing(epfl_word_dictionary, epfl_word_category)

print(epfl_word_pair)

# Applying the classify_text function to each text
epfl_text <- epfl_text %>%
  mutate(category = sapply(text, classify_text, word_pair = epfl_word_pair))

print(epfl_text)

goe

goe_text <- subset(data_posts_content, id == "goe.csv")
goe_corpus <- corpus_maker(goe_text)

goe_clean_text <- cleanup(goe_corpus, c("can", "new", "amp", "will"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
goe_word_dictionary <- freq_terms(goe_clean_text, 15)

goe_word_category <- c("Research", "Research", "Image", "Education", "Engagement", "Image", "Engagement", "Education")

goe_word_pair <- pairing(goe_word_dictionary, goe_word_category)

print(goe_word_pair)

# Applying the classify_text function to each text
goe_text <- goe_text %>%
  mutate(category = sapply(text, classify_text, word_pair = goe_word_pair))

print(goe_text)

harvard

harvard_text <- subset(data_posts_content, id == "harvard.csv")
harvard_corpus <- corpus_maker(harvard_text)

harvard_clean_text <- cleanup(harvard_corpus, c("new", "can", "will", "summer", "year", "first", "may", "-", "recent", "years", "one", "said", "time", "many", "world", "change"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
harvard_word_dictionary <- freq_terms(harvard_clean_text, 50)

harvard_word_category <- c("Image", "Education", "Research", "Image", "Research", "Research", "Engagement", "Education", "Society", "Education", "Society", "Research", "Education", "Image", "Research", NA, "Research", "Education", "Education", NA, "Research", "Society", "Society", "Research", "Education")

harvard_word_pair <- pairing(harvard_word_dictionary, harvard_word_category)

print(harvard_word_pair)

harvard_text <- harvard_text %>%
  mutate(category = sapply(text, classify_text, word_pair = harvard_word_pair))

print(harvard_text)

leicester

leicester_text <- subset(data_posts_content, id == "leicester.csv")
leicester_corpus <- corpus_maker(leicester_text)

leicester_clean_text <- cleanup(leicester_corpus, c("👉", "day", "new", "clear", "year", "can", "will", "time", "space", "one", "first"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
leicester_word_dictionary <- freq_terms(leicester_clean_text, 60)

leicester_word_category <- c("Image", NA, "Image", "Engagement", "Society", "Education", "Education", "Research", "Engagement", "Engagement", "Image", "Engagement", "Education", "Engamement", "Engagement", "Education", "Image", "Engagement", "Education", NA, "Image")

leicester_word_pair <- pairing(leicester_word_dictionary, leicester_word_category)

print(leicester_word_pair)

leicester_text <- leicester_text %>%
  mutate(category = sapply(text, classify_text, word_pair = leicester_word_pair))

print(leicester_text)

manchester

manchester_text <- subset(data_posts_content, id == "manchester.csv")
manchester_corpus <- corpus_maker(manchester_text)

manchester_clean_text <- cleanup(manchester_corpus, c("👇", "can", "just", "will", "get", "well", "one", "help", "now", "new", "read", "congratulations"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
manchester_word_dictionary <- freq_terms(manchester_clean_text, 50)

manchester_word_category <- c(NA, "Image", "Engagement", "Image", "Education", "Education", "Engagement", "Image", "Education", "Society", "Research", "Research", "Education", "Education", "Education", "Engagement", "Research", "Education", "Research")

manchester_word_pair <- pairing(manchester_word_dictionary, manchester_word_category)

print(manchester_word_pair)

manchester_text <- manchester_text %>%
  mutate(category = sapply(text, classify_text, word_pair = manchester_word_pair))

print(manchester_text)

mit

mit_text <- subset(data_posts_content, id == "mit.csv")
mit_corpus <- corpus_maker(mit_text)

mit_clean_text <- cleanup(mit_corpus, c("new", "can", "says", "“", "’", "—", "", "may", "first", "will", "using", "way", "one", "science"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
mit_word_dictionary <- freq_terms(mit_clean_text, 40)

mit_word_category <- c("Image", "Research", NA, "Image", NA, "Education", "Research", NA, "Education", "Research", "Research", "Education", "Research", "Research", "Image", "Education", "Education", "Research", "Research", "Research", "Research", "Research", "Research", "Research", "Research", "Society")

mit_word_pair <- pairing(mit_word_dictionary, mit_word_category)

print(mit_word_pair)

mit_text <- mit_text %>%
  mutate(category = sapply(text, classify_text, word_pair = mit_word_pair))

print(mit_text)

sb

sb_text <- subset(data_posts_content, id == "sb.csv")
sb_corpus <- corpus_maker(sb_text)

sb_clean_text <- cleanup(sb_corpus, c("new", "—", "will", "week", "can", "amp", "via", "now", "future", "first"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
sb_word_dictionary <- freq_terms(sb_clean_text, 30)

sb_word_category <- c("Image", "Research", "Image", "Image", "Image", NA, "Research", "Image", "Image", "Research", "Education", "Image", "Education", "Image", "Education", "Image", "Society")

sb_word_pair <- pairing(sb_word_dictionary, sb_word_category)

print(sb_word_pair)

sb_text <- sb_text %>%
  mutate(category = sapply(text, classify_text, word_pair = sb_word_pair))

print(sb_text)

stanford

stanford_text <- subset(data_posts_content, id == "stanford.csv")
stanford_corpus <- corpus_maker(stanford_text)

stanford_clean_text <- cleanup(stanford_corpus, c("new", "will"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
stanford_word_dictionary <- freq_terms(stanford_clean_text, 25)

stanford_word_category <- c("Image", "Image", "Research", "Education", NA)

stanford_word_pair <- pairing(stanford_word_dictionary, stanford_word_category)

print(stanford_word_pair)

stanford_text <- stanford_text %>%
  mutate(category = sapply(text, classify_text, word_pair = stanford_word_pair))

print(stanford_text)

trinity

trinity_text <- subset(data_posts_content, id == "trinity.csv")
trinity_corpus <- corpus_maker(trinity_text)

trinity_clean_text <- cleanup(trinity_corpus, c("amp", "read", "new", "can", "will", "week", "work", "great", "day", "visit", "irish", "first", "led", "congratulations"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
trinity_word_dictionary <- freq_terms(trinity_clean_text, 40)
Warning: empty document(s): 754
trinity_word_category <- c("Image", "Education", "Research", "Image", "Research", "Image", "Research", "Society", "Society", "Education", "Engagement", "Engagement", "Education", "Image", "Education", "Image", "Research")

trinity_word_pair <- pairing(trinity_word_dictionary, trinity_word_category)

print(trinity_word_pair)

trinity_text <- trinity_text %>%
  mutate(category = sapply(text, classify_text, word_pair = trinity_word_pair))

print(trinity_text)

wv

wv_text <- subset(data_posts_content, id == "wv.csv")
wv_corpus <- corpus_maker(wv_text)

wv_clean_text <- cleanup(wv_corpus, c("💛💙", "🙌", "👉", "happy", "day", "see", "week", "great", "can", "will", "well", "know", "new", "now", "get", "just"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
wv_word_dictionary <- freq_terms(wv_clean_text, 40)

wv_word_category <- c(NA, "Image", "Image", NA, "Image", "Image", "Engagement", "Education", "Education", "Education", NA)

wv_word_pair <- pairing(wv_word_dictionary, wv_word_category)

print(wv_word_pair)

wv_text <- wv_text %>%
  mutate(category = sapply(text, classify_text, word_pair = wv_word_pair))

print(wv_text)

yale

yale_text <- subset(data_posts_content, id == "yale.csv")
yale_corpus <- corpus_maker(yale_text)

yale_clean_text <- cleanup(yale_corpus, c("new", "—", "will", "can", "'", "first", "work", "read", "help", "year"))
Warning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documentsWarning: transformation drops documents
yale_word_dictionary <- freq_terms(yale_clean_text, 70)

yale_word_category <- c("Research", "Education", "Education", NA, "Research", "Research", "Research", "Research", "Education", "Image", "Research", "Education", "Image", "Research", NA, "Research", "Image", "Image", "Research", "Image", "Research", "Society", "Research", "Society", "Society")

yale_word_pair <- pairing(yale_word_dictionary, yale_word_category)

print(yale_word_pair)

yale_text <- yale_text %>%
  mutate(category = sapply(text, classify_text, word_pair = yale_word_pair))

print(yale_text)

Sentiment Analysis

Function for emotions

emotions_maker <- function(texts, hei_name){
  text_emotion <- get_nrc_sentiment(texts)
  
  # Proportions for text values 
  emotion_proportions <- colSums(prop.table(text_emotion[, 1:8]))
  
  # Visualization of percentages on each emotion found on posts
  barplot <- barplot(
    sort(colSums(prop.table(text_emotion[, 1:8]))), 
    horiz = TRUE, 
    cex.names = 0.7, 
    las = 1, 
    main = sprintf("Emotions found in %s's texts", hei_name), 
    xlab="Percentage",
    xlim = c(0, max(emotion_proportions) * 1.2)
  )
  
  text(
    x = sort(emotion_proportions),
    y = barplot,                    
    labels = sprintf("%.2f%%", 100 * sort(emotion_proportions)),
    pos = 4,                               
    cex = 0.7                              
  )
}

Duke

duke_text$sentiment <- round(get_sentiment(duke_text$text), 2)

print(duke_text)
emotions_maker(duke_text$text, "duke")

epfl

epfl_text$sentiment <- round(get_sentiment(epfl_text$text), 2)

print(epfl_text)
emotions_maker(epfl_text$text, "epfl")

goe

goe_text$sentiment <- round(get_sentiment(goe_text$text), 2)

print(goe_text)
emotions_maker(goe_text$text, "goe")

harvard

harvard_text$sentiment <- round(get_sentiment(harvard_text$text), 2)

print(harvard_text)
emotions_maker(harvard_text$text, "harvard")

leicester

leicester_text$sentiment <- round(get_sentiment(leicester_text$text), 2)

print(leicester_text)
emotions_maker(leicester_text$text, "leicester")

manchester

manchester_text$sentiment <- round(get_sentiment(manchester_text$text), 2)

print(manchester_text)
emotions_maker(manchester_text$text, "manchester")

mit

mit_text$sentiment <- round(get_sentiment(mit_text$text), 2)

print(mit_text)
emotions_maker(mit_text$text, "mit")

sb

sb_text$sentiment <- round(get_sentiment(sb_text$text), 2)

print(sb_text)
emotions_maker(sb_text$text, "sb")

stanford

stanford_text$sentiment <- round(get_sentiment(stanford_text$text), 2)

print(stanford_text)
emotions_maker(stanford_text$text, "stanford")

trinity

trinity_text$sentiment <- round(get_sentiment(trinity_text$text), 2)

print(trinity_text)
emotions_maker(trinity_text$text, "trinity")

wv

wv_text$sentiment <- round(get_sentiment(wv_text$text), 2)

print(wv_text)
emotions_maker(wv_text$text, "wv")

yale

yale_text$sentiment <- round(get_sentiment(yale_text$text), 2)

print(yale_text)
emotions_maker(yale_text$text, "yale")

Creation of average sentiment for each HEI

# Combine all the sentiment scores into one data frame
all_sentiments <- rbind(
    duke_text[, c("id", "sentiment")],
    epfl_text[, c("id", "sentiment")],
    goe_text[, c("id", "sentiment")],
    harvard_text[, c("id", "sentiment")],
    leicester_text[, c("id", "sentiment")],
    manchester_text[, c("id", "sentiment")],
    mit_text[, c("id", "sentiment")],
    sb_text[, c("id", "sentiment")],
    stanford_text[, c("id", "sentiment")],
    trinity_text[, c("id", "sentiment")],
    wv_text[, c("id", "sentiment")],
    yale_text[, c("id", "sentiment")]
)

# Calculate the average sentiment for each HEI
hei_average_sentiments <- aggregate(sentiment ~ id, data = all_sentiments, FUN = mean) %>%
    rename(average_sentiment = sentiment)

# Print the table
print(hei_average_sentiments)

Clusters

Function for cosine matrix

cosine_matrix_maker <- function(table){
  numerical_columns <- table %>%
    select(-id)
  
  # Normalizing columns
  normalized_columns <- as.data.frame(scale(numerical_columns))
  
  # Transposing the data to compute similarity between columns
  transposed_data <- t(normalized_columns)
  
  # Computing the cosine similarity matrix
  similarity_matrix <- as.matrix(proxy::dist(transposed_data, method = "cosine"))
  
  # Converting distance to similarity
  similarity_matrix <- 1 - similarity_matrix
  
  # Print the similarity matrix
  print(similarity_matrix)
}  

Creating table for cluster algorithms

# Joining attribute percentage_tweets (percentage of tweets out of all posts) and percentage_replies (percentage of replies out of all posts) from number_posts also adding unique_hashtags (number of unique hashtags) and hashtag_percentage (percentage of posts that contain a hashtag) from hashtags, per HEI
cluster_table <- merge(select(hashtags, id, unique_hashtags, hashtag_percentage), select(data_ratio, id, percentage_tweets, percentage_replies), by = "id", all=TRUE)

# Joining attribute avg_posts_per_days (average of posts per day) from posts_per_day per HEI
cluster_table <- merge(cluster_table, select(posts_per_day, id, avg_posts_per_days), by = "id", all=TRUE)

# Joining attribute avg_posts_per_weeks (average of posts per week) from posts_per_week per HEI
cluster_table <- merge(cluster_table, select(posts_per_week, id, avg_posts_per_weeks), by = "id", all=TRUE)

# Joining attribute avg_posts_in_academic_time (average of posts during academic time)  from data_posts_academic per HEI
cluster_table <- merge(cluster_table, select(data_posts_academic, id, avg_posts_in_academic_time), by = "id", all=TRUE)

# Joining attribute avg_posts_in_vacation_time (average of posts during vacation time) from data_posts_vacations per HEI
cluster_table <- merge(cluster_table, select(data_posts_vacations, id, avg_posts_in_vacation_time), by = "id", all=TRUE)

# Joining attribute time_of_day_value (numerical value referring to time of day where every HEI made more posts) from favourite_hour_hei per HEI
cluster_table <- merge(cluster_table, select(favourite_hour_hei, id, time_of_day_value), by = "id", all=TRUE)

# Joining attribute url_percentage (percentage of posts that contain an url) from url_usage per HEI
cluster_table <- merge(cluster_table, select(url_usage, id, url_percentage), by = "id", all=TRUE)

# Joining attribute average_num_words (average number of words in the posts) from data_posts_content_metrics per HEI
cluster_table <- merge(cluster_table, select(data_posts_content_metrics, id, average_num_words), by = "id", all=TRUE)

# Joining attribute average_sentiment (average sentiment of posts) from hei_average_sentiments per HEI
cluster_table <- merge(cluster_table, select(hei_average_sentiments, id, average_sentiment), by = "id", all=TRUE)

print(cluster_table)
cosine_matrix_maker(cluster_table)
                           unique_hashtags hashtag_percentage percentage_tweets percentage_replies avg_posts_per_days avg_posts_per_weeks
unique_hashtags                 1.00000000        0.714423178        0.16552198        -0.15047782         0.09423180          0.09460315
hashtag_percentage              0.71442318        1.000000000        0.13470812        -0.34565768        -0.01048208         -0.01015282
percentage_tweets               0.16552198        0.134708118        1.00000000        -0.73554060         0.30916774          0.30775916
percentage_replies             -0.15047782       -0.345657681       -0.73554060         1.00000000         0.09388665          0.09480553
avg_posts_per_days              0.09423180       -0.010482084        0.30916774         0.09388665         1.00000000          0.99999797
avg_posts_per_weeks             0.09460315       -0.010152816        0.30775916         0.09480553         0.99999797          1.00000000
avg_posts_in_academic_time      0.12447464       -0.033819668        0.32273419         0.06234258         0.97881722          0.97877880
avg_posts_in_vacation_time      0.09368691        0.004006656        0.16186441         0.20825601         0.97718481          0.97729844
time_of_day_value              -0.22814041       -0.175920221        0.30874190        -0.21552038        -0.05609489         -0.05648404
url_percentage                  0.09757172        0.395409742        0.52625841        -0.84342954         0.06010475          0.05920162
average_num_words               0.46150077        0.713419420        0.06039969        -0.48123350        -0.36976005         -0.36993229
average_sentiment               0.59299320        0.645971142       -0.31588590        -0.04093998        -0.49647197         -0.49605792
                           avg_posts_in_academic_time avg_posts_in_vacation_time time_of_day_value url_percentage average_num_words average_sentiment
unique_hashtags                            0.12447464                0.093686911       -0.22814041     0.09757172        0.46150077        0.59299320
hashtag_percentage                        -0.03381967                0.004006656       -0.17592022     0.39540974        0.71341942        0.64597114
percentage_tweets                          0.32273419                0.161864406        0.30874190     0.52625841        0.06039969       -0.31588590
percentage_replies                         0.06234258                0.208256011       -0.21552038    -0.84342954       -0.48123350       -0.04093998
avg_posts_per_days                         0.97881722                0.977184814       -0.05609489     0.06010475       -0.36976005       -0.49647197
avg_posts_per_weeks                        0.97877880                0.977298444       -0.05648404     0.05920162       -0.36993229       -0.49605792
avg_posts_in_academic_time                 1.00000000                0.953796153       -0.14950356     0.07310820       -0.37164829       -0.49875096
avg_posts_in_vacation_time                 0.95379615                1.000000000       -0.07107087     0.02168275       -0.32110128       -0.41192255
time_of_day_value                         -0.14950356               -0.071070873        1.00000000    -0.01251319       -0.10921047       -0.29227459
url_percentage                             0.07310820                0.021682748       -0.01251319     1.00000000        0.63011343        0.17847188
average_num_words                         -0.37164829               -0.321101278       -0.10921047     0.63011343        1.00000000        0.84649940
average_sentiment                         -0.49875096               -0.411922553       -0.29227459     0.17847188        0.84649940        1.00000000

Based on this cosine matrix we decided to remove average_num_words (high similarity to unique_hashtags, hashtag_percentage and url_percentage), unique_hashtags (high similarity to hashtag_percentage and the distances of hashtag_percentage to the other columns is lower) and avg_posts_per_weeks, avg_posts_in_academic_time, avg_posts_in_vacation_time are also all removed (due to their similarity to avg_posts_per_days)

cluster_table <- cluster_table %>%
  select(-average_num_words, -unique_hashtags, -avg_posts_per_weeks, -avg_posts_in_academic_time, -avg_posts_in_vacation_time)

cosine_matrix_maker(cluster_table)
                   hashtag_percentage percentage_tweets percentage_replies avg_posts_per_days time_of_day_value url_percentage average_sentiment
hashtag_percentage         1.00000000         0.1347081        -0.34565768        -0.01048208       -0.17592022     0.39540974        0.64597114
percentage_tweets          0.13470812         1.0000000        -0.73554060         0.30916774        0.30874190     0.52625841       -0.31588590
percentage_replies        -0.34565768        -0.7355406         1.00000000         0.09388665       -0.21552038    -0.84342954       -0.04093998
avg_posts_per_days        -0.01048208         0.3091677         0.09388665         1.00000000       -0.05609489     0.06010475       -0.49647197
time_of_day_value         -0.17592022         0.3087419        -0.21552038        -0.05609489        1.00000000    -0.01251319       -0.29227459
url_percentage             0.39540974         0.5262584        -0.84342954         0.06010475       -0.01251319     1.00000000        0.17847188
average_sentiment          0.64597114        -0.3158859        -0.04093998        -0.49647197       -0.29227459     0.17847188        1.00000000

Function for cluster method

cluster_maker <- function(num_clusters, table){
  # Excluding id column for clustering
  cluster_data <- select(table, -id)
    
  # Scaling the data for kmeans method
  scaled_data <- scale(cluster_data)
  
  kmeans_model <- kmeans(scaled_data, centers = num_clusters, nstart = 10)

  # Extract cluster assignments
  cluster_assignments <- kmeans_model$cluster
  
  # Create a data frame combining original data with cluster assignments
  clustered_data <- cbind(cluster_table$id, cluster_data, cluster = cluster_assignments)
  
  clustered_data <- clustered_data[, c("cluster_table$id", "cluster")]
  
  print(clustered_data)
}

Function to discover best number of clusters

elbow_maker <- function(table){
  cluster_data <- select(table, -id)
  scaled_data <- scale(cluster_data)
  
  wss <- vector()
  range <- 1:10
  
  for (k in range) {
    kmeans_model <- kmeans(scaled_data, centers = k, nstart = 10)
    wss[k] <- kmeans_model$tot.withinss
  }
  
  elbow_df <- data.frame(k = range, WSS = wss)
  ggplot(elbow_df, aes(x = k, y = WSS)) +
    geom_line() +
    geom_point() +
    labs(x = "Number of Clusters", y = "Within-Cluster Sum of Squares (WCSS)",
         title = "Elbow Method for Optimal k") +
    theme_minimal()
}

Plot of Elbow Method and selection of best number of cluster to view how HEIs are grouped

elbow_maker(cluster_table)

cluster_maker(4, cluster_table)
LS0tCnRpdGxlOiAiUHJvamVjdCIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQojIyMgRm9yIFIgYmVnaW5uZXJzCk5ldyBjaHVuayAqQ3RybCtBbHQrSSoKCkV4ZWN1dGUgY2h1bmsgKkN0cmwrU2hpZnQrRW50ZXIqCgpFeGVjdXRlIGFsbCBjaHVua3MgKkN0cmwrQWx0K1IqCgpIVE1MIHByZXZpZXcgKkN0cmwrU2hpZnQrSyoKCiMgTGlicmFyeSBwcmVwYXJhdGlvbnMKCmBgYHtyfQpsaWJyYXJ5KHJlYWRyKQpsaWJyYXJ5KGRwbHlyKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKbGlicmFyeShnZ3Bsb3QyKQpsaWJyYXJ5KHJlc2hhcGUyKQpsaWJyYXJ5KHN0YXRzKQpsaWJyYXJ5KHRtKQpsaWJyYXJ5KHRleHQydmVjKQpsaWJyYXJ5KHRleHRzdGVtKQpsaWJyYXJ5KHN5dXpoZXQpCmBgYAoKIyBEYXRhIEltcG9ydAoKYGBge3J9CmRhdGFfcG9zdHMgPC0gcmVhZC5jc3YoIn4vNHllYXIvMnNlbWVzdGVyL2R0SUkvQ1NWcy9IRUlzLmNzdiIsICAgICAjIFBsYWNlIGhlcmUgeW91ciBwYXRoIHRvIGZpbGUKICAgICAgICAgICAgICAgICBjb2xDbGFzc2VzID0gYyh0d2VldF9pZCA9ICJjaGFyYWN0ZXIiKSkKCiMgTW9kaWZ5aW5nIGNyZWF0ZWRfYXQgdHlwZSBzbyB0aGF0IGF0dHJpYnV0ZSBjYW4gYmUgdXNlZCBtb3JlIGVhc2lseSAKZGF0YV9wb3N0cyRjcmVhdGVkX2F0IDwtIGFzLlBPU0lYY3QoZGF0YV9wb3N0cyRjcmVhdGVkX2F0LAogICAgICAgICAgICAgICAgICAgICAgICAgICAgICBmb3JtYXQ9ICIlWS0lbS0lZFQlSDolTTolUyIsIHR6PSJVVEMiKQoKI1ZpZXcoZGF0YSkKc3VtbWFyeShkYXRhX3Bvc3RzKQpgYGAKCiMgRGF0YSBQcmUtUHJvY2Vzc2luZwoKYGBge3J9CiMgQ291bnQgb2YgaG93IG1hbnkgZW50cmllcyBlYWNoIEhFSSBoYXMKbnVtYmVyX3Bvc3RzIDwtIGRhdGFfcG9zdHMgJT4lCiAgICAgICAgICAgICAgZ3JvdXBfYnkoaWQpICU+JSBzdW1tYXJpc2UoY291bnQgPSBuKCkpCgpudW1iZXJfcG9zdHMKYGBgCgojIyMgU2luY2UgY29tcGx1dGVuc2Ugb25seSBoYXMgMSBlbnRyeSB3ZSBjYW4ndCBsZWFybiBhbnl0aGluZyBmcm9tIGl0LCBzbyB3ZSByZW1vdmVkIGl0CgpgYGB7cn0KZGF0YV9wb3N0cyA8LSBkYXRhX3Bvc3RzW2RhdGFfcG9zdHMkaWQgIT0gImNvbXBsdXRlbnNlLmNzdiIsIF0KYGBgCgojIyMgVmlzdWFsaXphdGlvbiBvZiBhbGwgcG9zdHMsIGp1c3QgdHdlZXRzIGFuZCBqdXN0IHJlcGxpZXMKCmBgYHtyfQpudW1iZXJfcG9zdHMgPC0gZGF0YV9wb3N0cyAlPiUKICAgICAgICAgICAgICBncm91cF9ieShpZCkgJT4lIHN1bW1hcmlzZShwb3N0cyA9IG4oKSkKCm51bWJlcl90d2VldHMgPC0gZGF0YV9wb3N0c1tkYXRhX3Bvc3RzJHR5cGUgPT0gIlR3ZWV0IiwgXSAlPiUKICAgICAgICAgICAgICBncm91cF9ieShpZCkgJT4lIHN1bW1hcmlzZSh0d2VldHMgPSBuKCkpCgpudW1iZXJfcmVwbGllcyA8LSBkYXRhX3Bvc3RzW2RhdGFfcG9zdHMkdHlwZSA9PSAiUmVwbHkiLCBdICU+JQogICAgICAgICAgICAgIGdyb3VwX2J5KGlkKSAlPiUgc3VtbWFyaXNlKHJlcGxpZXMgPSBuKCkpCgpwcmludChudW1iZXJfcG9zdHMpCnByaW50KG51bWJlcl90d2VldHMpCnByaW50KG51bWJlcl9yZXBsaWVzKQpgYGAKCiMjIyBDYWxjdWxhdGluZyB0aGUgcGVyY2VudGFnZSBvZiB0d2VldHMgYW5kIHJlcGxpZXMgYmFzZWQgb24gYWxsIHBvc3RzCgpgYGB7cn0KIyBNZXJnaW5nIHRoZSBjb3VudHMgb2YgdHdlZXRzIGFuZCByZXBsaWVzIHdpdGggdGhlIGNvdW50IG9mIHBvc3RzCmRhdGFfcmF0aW8gPC0gbWVyZ2UobnVtYmVyX3Bvc3RzLCBudW1iZXJfdHdlZXRzLCBieSA9ICJpZCIsIGFsbCA9IFRSVUUpCmRhdGFfcmF0aW8gPC0gbWVyZ2UoZGF0YV9yYXRpbywgbnVtYmVyX3JlcGxpZXMsIGJ5ID0gImlkIiwgYWxsID0gVFJVRSkKCgpkYXRhX3JhdGlvJHBlcmNlbnRhZ2VfdHdlZXRzIDwtIHJvdW5kKCgoZGF0YV9yYXRpbyR0d2VldHMgLyBkYXRhX3JhdGlvJHBvc3RzKSAqIDEwMCksIDIpCmRhdGFfcmF0aW8kcGVyY2VudGFnZV9yZXBsaWVzIDwtIHJvdW5kKCgoZGF0YV9yYXRpbyRyZXBsaWVzIC8gZGF0YV9yYXRpbyRwb3N0cykgKiAxMDApLCAyKQoKZGF0YV9yYXRpbyA8LSBkYXRhX3JhdGlvWywgYygiaWQiLCAicGVyY2VudGFnZV90d2VldHMiLCAicGVyY2VudGFnZV9yZXBsaWVzIildCgpkYXRhX3JhdGlvJHBlcmNlbnRhZ2VfcmVwbGllc1tpcy5uYShkYXRhX3JhdGlvJHBlcmNlbnRhZ2VfcmVwbGllcyldIDwtIDAKCnByaW50KGRhdGFfcmF0aW8pCmBgYAoKIyMgTkEgcmVtb3ZhbAoKIyMjIEZ1bmN0aW9uIHRvIHZpc3VhbGl6ZSB0aGUgbnVtYmVyIG9mIE5BcyBpbiBhbGwgY29sdW1ucwoKYGBge3J9Cm5hX2NvdW50IDwtIGZ1bmN0aW9uKCl7CiAgIyBDb3VudGluZyB0aGUgbnVtYmVyIG9mIE5BIHZhbHVlcyBmb3IgZWFjaCBjb2x1bW4KICBuYV9jb3VudCA8LSBjb2xTdW1zKGlzLm5hKGRhdGFfcG9zdHMpKQogIAogICMgQ3JlYXRpbmcgYSBuZXcgZGF0YSBmcmFtZSB3aXRoIHRoZSBOQSBjb3VudHMKICBuYV9jb3VudHNfdGFibGUgPC0gZGF0YS5mcmFtZShDb2x1bW4gPSBuYW1lcyhuYV9jb3VudCksIE5BX0NvdW50ID0gbmFfY291bnQpCiAgCiAgcHJpbnQobmFfY291bnRzX3RhYmxlKQp9CmBgYAoKIyMjIENhbGN1bGF0aW9ucyBvZiB2aWV3LCBmYXZvdXJpdGUsIHJldHdlZXQgYW5kIHJlcGx5IHBlcmNlbnRpbGVzIGFuZCB2aXN1YWxpemF0aW9uIG9mIE5BcyBpbiBhbGwgY29sdW1ucwoKYGBge3J9CmRhdGFfcG9zdHMgPC0gZGF0YV9wb3N0cyAlPiUKICBncm91cF9ieShpZCkgJT4lCiAgbXV0YXRlKHZpZXdfcGVyY2VudGlsZSA9IG50aWxlKHZpZXdfY291bnQsIDEwMCksCiAgICAgICAgIGZhdm9yaXRlX3BlcmNlbnRpbGUgPSBudGlsZShmYXZvcml0ZV9jb3VudCwgMTAwKSwKICAgICAgICAgcmV0d2VldF9wZXJjZW50aWxlID0gbnRpbGUocmV0d2VldF9jb3VudCwgMTAwKSwKICAgICAgICAgcmVwbHlfcGVyY2VudGlsZSA9IG50aWxlKHJlcGx5X2NvdW50LCAxMDApKSAlPiUKICByb3d3aXNlKCkgJT4lCiAgbXV0YXRlKGF2Z19wZXJjZW50aWxlID0gcm91bmQobWVhbihjKHZpZXdfcGVyY2VudGlsZSwgZmF2b3JpdGVfcGVyY2VudGlsZSwgcmV0d2VldF9wZXJjZW50aWxlLCByZXBseV9wZXJjZW50aWxlKSwgbmEucm0gPSBUUlVFKSwgMikpCgpuYV9jb3VudCgpCgpkYXRhX3BlcmNlbnRpbGUgPC0gZGF0YV9wb3N0c1ssIGMoImlkIiwgInZpZXdfcGVyY2VudGlsZSIsICJmYXZvcml0ZV9wZXJjZW50aWxlIiwgInJldHdlZXRfcGVyY2VudGlsZSIsICJyZXBseV9wZXJjZW50aWxlIiwgImF2Z19wZXJjZW50aWxlIildCgpwcmludChkYXRhX3BlcmNlbnRpbGUpCmBgYAoKIyMjIENhbGN1bGF0aW9uIG9mIHRoZSBtYXhpbXVtIG51bWJlciBvZiB2aWV3cyBmb3IgZWFjaCBIRUkKCmBgYHtyfQptYXhfdmlld19jb3VudHMgPC0gdGFwcGx5KGRhdGFfcG9zdHMkdmlld19jb3VudCwgZGF0YV9wb3N0cyRpZCwgbWF4LCBuYS5ybSA9IFRSVUUpCgpwcmludChtYXhfdmlld19jb3VudHMpCmBgYAoKIyMjIFJlbW92YWwgb2YgTkFzCgpgYGB7cn0KIyBGcm9tIHZpZXcgY291bnQKZGF0YV9wb3N0cyR2aWV3X2NvdW50IDwtIGlmZWxzZSgKICBpcy5uYShkYXRhX3Bvc3RzJHZpZXdfY291bnQpLAogIHJvdW5kKG1heF92aWV3X2NvdW50c1tkYXRhX3Bvc3RzJGlkXSAqIChkYXRhX3Bvc3RzJGF2Z19wZXJjZW50aWxlIC8gMTAwKSksCiAgZGF0YV9wb3N0cyR2aWV3X2NvdW50KQoKIyBGcm9tIHZpZXcgcGVyY2VudGlsZQpkYXRhX3Bvc3RzJHZpZXdfcGVyY2VudGlsZSA8LSBpZmVsc2UoCiAgaXMubmEoZGF0YV9wb3N0cyR2aWV3X3BlcmNlbnRpbGUpLAogIGRhdGFfcG9zdHMkYXZnX3BlcmNlbnRpbGUsCiAgZGF0YV9wb3N0cyR2aWV3X3BlcmNlbnRpbGUpCmBgYAoKIyMjIFZpc3VhbGl6YXRpb24gb2YgTkFzIGluIGFsbCBjb2x1bW5zCgpgYGB7cn0KbmFfY291bnQoKQpgYGAKCiMgRGF0YSBQcm9jZXNzaW5nCgojIyMgRnVuY3Rpb24gdG8gY2FsY3VsYXRlIGF2ZXJhZ2UgcG9zdHMKCmBgYHtyfQphdmVyYWdlX3Bvc3RzIDwtIGZ1bmN0aW9uKHRpbWVmcmFtZSl7CiAgIyBDYWxjdWxhdGlvbiBvZiB0aGUgdGltZWZyYW1lIGJldHdlZW4gZWFybGllc3QgYW5kIGxhdGVzdCBwb3N0IGZvciBlYWNoIEhFSQogIGRhdGVfcmFuZ2UgPC0gZGF0YV9wb3N0cyAlPiUKICAgIGdyb3VwX2J5KGlkKSAlPiUKICAgIHN1bW1hcmlzZShtaW5fZGF0ZSA9IG1pbihjcmVhdGVkX2F0KSwKICAgICAgICAgICAgICBtYXhfZGF0ZSA9IG1heChjcmVhdGVkX2F0KSkgJT4lCiAgICBtdXRhdGUobnVtX2RheXMgPSBhcy5udW1lcmljKGRpZmZ0aW1lKG1heF9kYXRlLCBtaW5fZGF0ZSwgdW5pdHMgPSB0aW1lZnJhbWUpKSkKICAKICAjIE5hbWluZyB0aGUgY29sdW1uIHJlc3BlY3RpbmcgdGhlIHRpbWVmcmFtZQogIGNvbHVtbl9uYW1lIDwtIHBhc3RlMCgiYXZnX3Bvc3RzX3Blcl8iLCB0aW1lZnJhbWUpCiAgCiAgIyBDYWxjdWxhdGlvbiBvZiB0aGUgbnVtYmVyIG9mIHBvc3RzIHBlciBkYXkgZm9yIGVhY2ggSEVJCiAgcG9zdHNfcGVyX3RpbWVmcmFtZSA8LSBudW1iZXJfcG9zdHMgJT4lCiAgICBsZWZ0X2pvaW4oZGF0ZV9yYW5nZSwgYnkgPSAiaWQiKSAlPiUKICAgIG11dGF0ZSghIWNvbHVtbl9uYW1lIDo9IHJvdW5kKChwb3N0cyAvIG51bV9kYXlzKSwgMikpCiAgCiAgcHJpbnQocG9zdHNfcGVyX3RpbWVmcmFtZSkKICByZXR1cm4ocG9zdHNfcGVyX3RpbWVmcmFtZSkKfQpgYGAKCmBgYHtyfQpwb3N0c19wZXJfZGF5IDwtIGF2ZXJhZ2VfcG9zdHMoImRheXMiKQpwb3N0c19wZXJfd2VlayA8LSBhdmVyYWdlX3Bvc3RzKCJ3ZWVrcyIpCmBgYAoKIyMjIFBsb3QgZm9yIHRoZSBhdmVyYWdlIG51bWJlciBvZiBwb3N0cyBwZXIgZGF5IGZvciBlYWNoIEhFSQoKYGBge3J9CmJhcnBsb3QocG9zdHNfcGVyX2RheSRhdmdfcG9zdHNfcGVyX2RheXMsCiAgICAgICAgbmFtZXMuYXJnID0gcG9zdHNfcGVyX2RheSRpZCwKICAgICAgICBtYWluID0gIkF2ZXJhZ2UgUG9zdHMgcGVyIERheSIsCiAgICAgICAgeGxhYiA9ICJIRUkiLAogICAgICAgIHlsYWIgPSAiQXZlcmFnZSBOdW1iZXIgb2YgUG9zdHMiLAogICAgICAgIHlsaW0gPSBjKDAsIG1heChwb3N0c19wZXJfZGF5JGF2Z19wb3N0c19wZXJfZGF5cykgKyAxKSwKICAgICAgICBsYXMgPSAyLAogICAgICAgIGNvbCA9ICIjMzQ5OERCIikKCiMgQWRkaW5nIHRleHQgbGFiZWxzIG92ZXIgZWFjaCBiYXIgYW5kIGFsaWduaW5nIGl0IHdpdGggdGhlIGNlbnRlciBvZiBlYWNoIGJhciAKdGV4dCh4ID0gYmFycGxvdChwb3N0c19wZXJfZGF5JGF2Z19wb3N0c19wZXJfZGF5cywgcGxvdCA9IEZBTFNFKSwKICAgICB5ID0gcG9zdHNfcGVyX2RheSRhdmdfcG9zdHNfcGVyX2RheXMsCiAgICAgbGFiZWxzID0gcm91bmQocG9zdHNfcGVyX2RheSRhdmdfcG9zdHNfcGVyX2RheXMsIDIpLAogICAgIHBvcyA9IDMpCmBgYAoKIyMjIFBsb3QgZm9yIHRoZSBhdmVyYWdlIG51bWJlciBvZiBwb3N0cyBwZXIgd2VlayBmb3IgZWFjaCBIRUkKCmBgYHtyfQpiYXJwbG90KHBvc3RzX3Blcl93ZWVrJGF2Z19wb3N0c19wZXJfd2Vla3MsCiAgICAgICAgbmFtZXMuYXJnID0gcG9zdHNfcGVyX3dlZWskaWQsCiAgICAgICAgbWFpbiA9ICJBdmVyYWdlIFBvc3RzIHBlciBXZWVrIiwKICAgICAgICB4bGFiID0gIkhFSSIsCiAgICAgICAgeWxhYiA9ICJBdmVyYWdlIE51bWJlciBvZiBQb3N0cyIsCiAgICAgICAgeWxpbSA9IGMoMCwgbWF4KHBvc3RzX3Blcl93ZWVrJGF2Z19wb3N0c19wZXJfd2Vla3MpICsgNSksCiAgICAgICAgbGFzID0gMiwKICAgICAgICBjb2wgPSAiI0U3NEMzQyIpCgp0ZXh0KHggPSBiYXJwbG90KHBvc3RzX3Blcl93ZWVrJGF2Z19wb3N0c19wZXJfd2Vla3MsIHBsb3QgPSBGQUxTRSksCiAgICAgeSA9IHBvc3RzX3Blcl93ZWVrJGF2Z19wb3N0c19wZXJfd2Vla3MsCiAgICAgbGFiZWxzID0gcm91bmQocG9zdHNfcGVyX3dlZWskYXZnX3Bvc3RzX3Blcl93ZWVrcywgMiksCiAgICAgcG9zID0gMykKYGBgCgojIyMgRGVmaW5pbmcgdGhlIGludGVydmFscyBvZiB0aW1lIGZvciB0aGUgYWNhZGVtaWMgeWVhcgoKYGBge3J9CmludGVydmFscyA8LSBsaXN0KAogIGludGVydmFsMSA9IGFzLlBPU0lYY3QoYygiMjAyMi0wOC0zMSIsICIyMDIyLTEyLTE1IikpLAogIGludGVydmFsMiA9IGFzLlBPU0lYY3QoYygiMjAyMy0wMS0wNCIsICIyMDIzLTA0LTAxIikpLAogIGludGVydmFsMyA9IGFzLlBPU0lYY3QoYygiMjAyMy0wNC0xNCIsICIyMDIzLTA2LTE1IikpCikKYGBgCgojIyMgRnVuY3Rpb24gdG8gY2hlY2sgaWYgYSBkYXRlIGZhbGxzIHdpdGhpbiBhIGdpdmVuIGludGVydmFsIG9mIHRpbWUgYW5kIGFwcGx5IGFwcHJvcHJpYXRlIEJvb2xlYW4KCmBgYHtyfQpjaGVja19pbnRlcnZhbCA8LSBmdW5jdGlvbihkYXRlKSB7CiAgZm9yIChpIGluIDE6bGVuZ3RoKGludGVydmFscykpIHsKICAgIGludGVydmFsX3N0YXJ0IDwtIGludGVydmFsc1tbaV1dWzFdCiAgICBpbnRlcnZhbF9lbmQgPC0gaW50ZXJ2YWxzW1tpXV1bMl0KICAgIGlmIChkYXRlID49IGludGVydmFsX3N0YXJ0ICYgZGF0ZSA8PSBpbnRlcnZhbF9lbmQpIHsKICAgICAgcmV0dXJuKFRSVUUpCiAgICB9CiAgfQogIHJldHVybihGQUxTRSkKfQpgYGAKCmBgYHtyfQpkYXRhX3Bvc3RzJGFjYWRlbWljX3llYXIgPC0gc2FwcGx5KGRhdGFfcG9zdHMkY3JlYXRlZF9hdCwgY2hlY2tfaW50ZXJ2YWwpCnByaW50KGRhdGEuZnJhbWUoaWQgPSBkYXRhX3Bvc3RzJGlkLCBhY2FkZW1pY195ZWFyID0gZGF0YV9wb3N0cyRhY2FkZW1pY195ZWFyKSkKYGBgCgojIyMgRnVuY3Rpb24gdG8gY291bnQgbnVtYmVyIG9mIHBvc3RzIGFuZCBhdmVyYWdlIHBlciBkYXkgZHVyaW5nIGFjYWRlbWljIHRpbWUgYW5kIHZhY2F0aW9uIHRpbWUKCmBgYHtyfQphbmFseXplX3Bvc3RzIDwtIGZ1bmN0aW9uKGFjYWRlbWljX3llYXJfZmlsdGVyKSB7CiAgIyBGaWx0ZXJpbmcgdGhlIGRhdGEgYmFzZWQgb24gdGhlIGFjYWRlbWljX3llYXJfZmlsdGVyCiAgZmlsdGVyZWRfZGF0YSA8LSBkYXRhX3Bvc3RzICU+JQogICAgZmlsdGVyKGFjYWRlbWljX3llYXIgPT0gYWNhZGVtaWNfeWVhcl9maWx0ZXIpCiAgCiAgIyBDb3VudCBvZiBkYXlzIGZvciBlYWNoIEhFSQogIHVuaXF1ZV9kYXlzIDwtIGZpbHRlcmVkX2RhdGEgJT4lCiAgICBncm91cF9ieShpZCkgJT4lCiAgICBzdW1tYXJpc2UodW5pcXVlX2RheXMgPSBuX2Rpc3RpbmN0KGFzLkRhdGUoY3JlYXRlZF9hdCkpKQogIAogICMgQ291bnQgb2YgcG9zdHMgZm9yIGVhY2ggSEVJCiAgbnVtYmVyX3Bvc3RzX2Jvb2xlYW4gPC0gZmlsdGVyZWRfZGF0YSAlPiUKICAgIGdyb3VwX2J5KGlkKSAlPiUKICAgIHN1bW1hcmlzZShjb3VudCA9IG4oKSkKICAKICAjIE5hbWluZyB0aGUgY29sdW1uIHJlc3BlY3RpbmcgdGhlIHRpbWUgcGVyaW9kCiAgdGltZSA8LSBpZmVsc2UoYWNhZGVtaWNfeWVhcl9maWx0ZXIsICJhY2FkZW1pY190aW1lIiwgInZhY2F0aW9uX3RpbWUiKQogIGNvbHVtbl9uYW1lIDwtIHBhc3RlMCgiYXZnX3Bvc3RzX2luXyIsIHRpbWUpCiAgCiAgIyBDb21iaW5hdGlvbiBvZiBkYXRhIGFuZCBjYWxjdWxhdGlvbiBvZiBhdmVyYWdlIHBvc3RzIHBlciBkYXkKICBjb21iaW5lZF9kYXRhIDwtIGxlZnRfam9pbih1bmlxdWVfZGF5cywgbnVtYmVyX3Bvc3RzX2Jvb2xlYW4sIGJ5ID0gImlkIikKICBjb21iaW5lZF9kYXRhIDwtIGNvbWJpbmVkX2RhdGEgJT4lCiAgICBtdXRhdGUoISFjb2x1bW5fbmFtZSA6PSByb3VuZCgoY291bnQgLyB1bmlxdWVfZGF5cyksIDIpKQogIAogIHByaW50KGNvbWJpbmVkX2RhdGEpCiAgcmV0dXJuKGNvbWJpbmVkX2RhdGEpCn0KYGBgCgpgYGB7cn0KZGF0YV9wb3N0c19hY2FkZW1pYyA8LSBhbmFseXplX3Bvc3RzKFRSVUUpCmRhdGFfcG9zdHNfdmFjYXRpb25zIDwtIGFuYWx5emVfcG9zdHMoRkFMU0UpCmBgYAoKIyMjIFBsb3QgZm9yIHRoZSBhdmVyYWdlIG51bWJlciBvZiBwb3N0cyBkdXJpbmcgYWNhZGVtaWMgdGltZSBmb3IgZWFjaCBIRUkKCmBgYHtyfQpiYXJwbG90KGRhdGFfcG9zdHNfYWNhZGVtaWMkYXZnX3Bvc3RzX2luX2FjYWRlbWljX3RpbWUsCiAgICAgICAgbmFtZXMuYXJnID0gZGF0YV9wb3N0c19hY2FkZW1pYyRpZCwKICAgICAgICBtYWluID0gIkF2ZXJhZ2UgUG9zdHMgZHVyaW5nIEFjYWRlbWljIFRpbWUiLAogICAgICAgIHhsYWIgPSAiSEVJIiwKICAgICAgICB5bGFiID0gIkF2ZXJhZ2UgTnVtYmVyIG9mIFBvc3RzIiwKICAgICAgICB5bGltID0gYygwLCBtYXgoZGF0YV9wb3N0c19hY2FkZW1pYyRhdmdfcG9zdHNfaW5fYWNhZGVtaWNfdGltZSkgKyA1KSwKICAgICAgICBsYXMgPSAyLAogICAgICAgIGNvbCA9ICIjMzQ0OTVFIikKCnRleHQoeCA9IGJhcnBsb3QoZGF0YV9wb3N0c19hY2FkZW1pYyRhdmdfcG9zdHNfaW5fYWNhZGVtaWNfdGltZSwgcGxvdCA9IEZBTFNFKSwKICAgICB5ID0gZGF0YV9wb3N0c19hY2FkZW1pYyRhdmdfcG9zdHNfaW5fYWNhZGVtaWNfdGltZSwKICAgICBsYWJlbHMgPSByb3VuZChkYXRhX3Bvc3RzX2FjYWRlbWljJGF2Z19wb3N0c19pbl9hY2FkZW1pY190aW1lLCAyKSwKICAgICBwb3MgPSAzKQpgYGAKCiMjIyBQbG90IGZvciB0aGUgYXZlcmFnZSBudW1iZXIgb2YgcG9zdHMgZHVyaW5nIHZhY2F0aW9uIHRpbWUgZm9yIGVhY2ggSEVJCgpgYGB7cn0KYmFycGxvdChkYXRhX3Bvc3RzX3ZhY2F0aW9ucyRhdmdfcG9zdHNfaW5fdmFjYXRpb25fdGltZSwKICAgICAgICBuYW1lcy5hcmcgPSBkYXRhX3Bvc3RzX3ZhY2F0aW9ucyRpZCwKICAgICAgICBtYWluID0gIkF2ZXJhZ2UgUG9zdHMgZHVyaW5nIFZhY2F0aW9uIFRpbWUiLAogICAgICAgIHhsYWIgPSAiSEVJIiwKICAgICAgICB5bGFiID0gIkF2ZXJhZ2UgTnVtYmVyIG9mIFBvc3RzIiwKICAgICAgICB5bGltID0gYygwLCBtYXgoZGF0YV9wb3N0c192YWNhdGlvbnMkYXZnX3Bvc3RzX2luX3ZhY2F0aW9uX3RpbWUpICsgNSksCiAgICAgICAgbGFzID0gMiwKICAgICAgICBjb2wgPSAiI0QzNTQwMCIpCgp0ZXh0KHggPSBiYXJwbG90KGRhdGFfcG9zdHNfdmFjYXRpb25zJGF2Z19wb3N0c19pbl92YWNhdGlvbl90aW1lLCBwbG90ID0gRkFMU0UpLAogICAgIHkgPSBkYXRhX3Bvc3RzX3ZhY2F0aW9ucyRhdmdfcG9zdHNfaW5fdmFjYXRpb25fdGltZSwKICAgICBsYWJlbHMgPSByb3VuZChkYXRhX3Bvc3RzX3ZhY2F0aW9ucyRhdmdfcG9zdHNfaW5fdmFjYXRpb25fdGltZSwgMiksCiAgICAgcG9zID0gMykKYGBgCgojIyBEYXRhIHByZXBhcmF0aW9uIGZvciBkYXRlcyAKCmBgYHtyfQojIENyZWF0aW5nIG5ldyB0YWJsZSB0aGF0IGNvbnRhaW5zIGEgbmV3IGNvbHVtbiBmb3IgdGhlIGRheSBvZiB0aGUgd2VlawpkYXRhX3Bvc3RzX2RheXMgPC0gZGF0YV9wb3N0cyAlPiUKICBtdXRhdGUoZGF5X29mX3dlZWsgPSB3ZWVrZGF5cyhjcmVhdGVkX2F0KSkKCiMgU2VsZWN0aW5nIG9ubHkgdGhlIGlkLCBjcmVhdGVkX2F0LCBhbmQgZGF5X29mX3dlZWsgY29sdW1ucyBmb3IgdGhlIG5ldyB0YWJsZQpkYXRhX3Bvc3RzX2RheXMgPC0gZGF0YV9wb3N0c19kYXlzICU+JQogIHNlbGVjdChpZCwgY3JlYXRlZF9hdCwgZGF5X29mX3dlZWspCgojIENyZWF0ZSBjb2x1bW4gaG91ciBmcm9tIGNyZWF0ZWRfYXQKZGF0YV9wb3N0c19kYXlzJGNyZWF0ZWRfaG91ciA8LSBhcy5udW1lcmljKGZvcm1hdChkYXRhX3Bvc3RzX2RheXMkY3JlYXRlZF9hdCwgIiVIIikpCgpwcmludChkYXRhX3Bvc3RzX2RheXMpCmBgYAoKYGBge3J9CiMgR3JvdXBpbmcgYnkgaWQgYW5kIGRheV9vZl93ZWVrLCB0aGVuIGNvdW50aW5nIHRoZSBudW1iZXIgb2YgcG9zdHMKbnVtYmVyX3Bvc3RzX2RheXMgPC0gZGF0YV9wb3N0c19kYXlzICU+JQogIGdyb3VwX2J5KGlkLCBkYXlfb2Zfd2VlaykgJT4lCiAgc3VtbWFyaXNlKGNvdW50ID0gbigpKQoKIyBHcm91cGluZyBieSBpZCwgZGF5X29mX3dlZWsgYW5kIGRheSBjcmVhdGVkIGF0LCB0aGVuIGNvdW50aW5nIHRoIGVudW1iZXIgb2YgdHdlZXRzCm51bWJlcl9wb3N0c19wZXJfZGF5IDwtIGRhdGFfcG9zdHNfZGF5cyAlPiUKICAgIG11dGF0ZShjcmVhdGVkX2RhdGUgPSBhcy5EYXRlKGNyZWF0ZWRfYXQpKSAlPiUgCiAgICBncm91cF9ieShpZCwgZGF5X29mX3dlZWssIGNyZWF0ZWRfZGF0ZSkgJT4lCiAgICBzdW1tYXJpemUoY291bnQgPSBuKCkpCgojIEZpbmRpbmcgZm9yIGVhY2ggSEVJIHRoZSBhdmVyYWdlIGNvdW50IG9mIHBvc3RzIHBlciBkYXkKYXZlcmFnZV9udW1iZXJfcG9zdHNfcGVyX2RheSA8LSBudW1iZXJfcG9zdHNfcGVyX2RheSAlPiUKICBncm91cF9ieShpZCwgZGF5X29mX3dlZWspICU+JQogIHN1bW1hcmlzZShhdmVyYWdlX2NvdW50ID0gcm91bmQobWVhbihjb3VudCksIDIpKQoKcHJpbnQobnVtYmVyX3Bvc3RzX2RheXMpCmBgYAoKIyMgSGlnaGVzdCBhbmQgbG93ZXN0IHBvc3RzCgpgYGB7cn0KIyBGaW5kaW5nIHRoZSBIRUkgd2l0aCB0aGUgbG93ZXN0IGNvdW50IG9mIHBvc3RzIHBlciBkYXkKbG93ZXN0X2NvdW50IDwtIG51bWJlcl9wb3N0c19kYXlzICU+JQogIGdyb3VwX2J5KGRheV9vZl93ZWVrKSAlPiUKICBzbGljZV9taW4ob3JkZXJfYnkgPSBjb3VudCkgJT4lCiAgc2VsZWN0KGRheV9vZl93ZWVrLCBpZCwgY291bnQpCgojIEZpbmRpbmcgdGhlIEhFSSB3aXRoIHRoZSBoaWdoZXN0IGNvdW50IG9mIHBvc3RzIHBlciBkYXkKaGlnaGVzdF9jb3VudCA8LSBudW1iZXJfcG9zdHNfZGF5cyAlPiUKICBncm91cF9ieShkYXlfb2Zfd2VlaykgJT4lCiAgc2xpY2VfbWF4KG9yZGVyX2J5ID0gY291bnQpICU+JQogIHNlbGVjdChkYXlfb2Zfd2VlaywgaWQsIGNvdW50KQoKIyBDb21iaW5lIHRoZSByZXN1bHRzCmhpZ2hfbG93X0hFSSA8LSBiaW5kX3Jvd3MobG93ZXN0X2NvdW50LCBoaWdoZXN0X2NvdW50KSAlPiUKICBhcnJhbmdlKGRheV9vZl93ZWVrKQoKcHJpbnQoaGlnaF9sb3dfSEVJKQpgYGAKCiMjIyBQbG90IGZvciB0aGUgaGlnaGVzdCBhbmQgbG93ZXN0IGNvdW50IG9mIHBvc3RzIHBlciBkYXkgZm9yIGVhY2ggZGF5IG9mIHRoZSB3ZWVrCgpgYGB7cn0KZ2dwbG90KGhpZ2hfbG93X0hFSSwgYWVzKHggPSBkYXlfb2Zfd2VlaywgeSA9IGNvdW50LCBmaWxsID0gaWQpKSArCiAgZ2VvbV9iYXIoc3RhdCA9ICJpZGVudGl0eSIsIHBvc2l0aW9uID0gImRvZGdlIikgKwogIGdlb21fdGV4dChhZXMobGFiZWwgPSBjb3VudCksCiAgICAgICAgICAgIHBvc2l0aW9uID0gcG9zaXRpb25fZG9kZ2Uod2lkdGggPSAwLjkpLAogICAgICAgICAgICB2anVzdCA9IC0wLjUsCiAgICAgICAgICAgIHNpemUgPSAzKSArCiAgbGFicyh0aXRsZSA9ICJIaWdoZXN0IGFuZCBMb3dlc3QgQ291bnQgb2YgUG9zdHMgcGVyIERheSBmb3IgRWFjaCBEYXkgb2YgdGhlIFdlZWsiLAogICAgICAgeCA9ICJEYXkgb2YgdGhlIFdlZWsiLCB5ID0gIkNvdW50IikgKwogIHNjYWxlX2ZpbGxfbWFudWFsKHZhbHVlcyA9IHJhaW5ib3cobGVuZ3RoKHVuaXF1ZShoaWdoX2xvd19IRUkkaWQpKSkpICsKICB0aGVtZV9taW5pbWFsKCkgKwogIHRoZW1lKGxlZ2VuZC50aXRsZSA9IGVsZW1lbnRfYmxhbmsoKSkKYGBgCgojIyBBdmVyYWdlIG9mIHBvc3RzCgpgYGB7cn0KIyBGaW5kaW5nIHRoZSBIRUkgd2l0aCBsb3dlc3QgYW5kIGhpZ2hlc3QgYXZlcmFnZWQgY291bnQgb2YgcG9zdHMgcGVyIGRheQpoaWdoX2xvd19hdmVyYWdlX0hFSXMgPC0gYXZlcmFnZV9udW1iZXJfcG9zdHNfcGVyX2RheSAlPiUKICBncm91cF9ieShkYXlfb2Zfd2VlaykgJT4lCiAgZmlsdGVyKGF2ZXJhZ2VfY291bnQgPT0gbWF4KGF2ZXJhZ2VfY291bnQpIHwgYXZlcmFnZV9jb3VudCA9PSBtaW4oYXZlcmFnZV9jb3VudCkpICU+JQogIGFycmFuZ2UoZGF5X29mX3dlZWssIGlmZWxzZShhdmVyYWdlX2NvdW50ID09IG1pbihhdmVyYWdlX2NvdW50KSwgYXZlcmFnZV9jb3VudCwgLWF2ZXJhZ2VfY291bnQpKQoKcHJpbnQoaGlnaF9sb3dfYXZlcmFnZV9IRUlzKQpgYGAKCiMjIyBQbG90IGZvciB0aGUgaGlnaGVzdCBhbmQgbG93ZXN0IGF2ZXJhZ2UgY291bnQgb2YgcG9zdHMgcGVyIGRheSBmb3IgZWFjaCBkYXkgb2YgdGhlIHdlZWsKCmBgYHtyfQpnZ3Bsb3QoaGlnaF9sb3dfYXZlcmFnZV9IRUlzLCBhZXMoeCA9IGRheV9vZl93ZWVrLCB5ID0gYXZlcmFnZV9jb3VudCwgZmlsbCA9IGlkKSkgKwogIGdlb21fYmFyKHN0YXQgPSAiaWRlbnRpdHkiLCBwb3NpdGlvbiA9ICJkb2RnZSIpICsKICBnZW9tX3RleHQoYWVzKGxhYmVsID0gcm91bmQoYXZlcmFnZV9jb3VudCwgMikpLAogICAgICAgICAgICBwb3NpdGlvbiA9IHBvc2l0aW9uX2RvZGdlKHdpZHRoID0gMC43KSwKICAgICAgICAgICAgdmp1c3QgPSAtMC41LAogICAgICAgICAgICBzaXplID0gMykgKwogIGxhYnModGl0bGUgPSAiSGlnaGVzdCBhbmQgTG93ZXN0IEF2ZXJhZ2UgQ291bnQgb2YgUG9zdHMgcGVyIERheSBmb3IgRWFjaCBEYXkgb2YgdGhlIFdlZWsiLAogICAgICAgeCA9ICJEYXkgb2YgdGhlIFdlZWsiLCB5ID0gIkF2ZXJhZ2UgQ291bnQiKSArCiAgc2NhbGVfZmlsbF9tYW51YWwodmFsdWVzID0gcmFpbmJvdyhsZW5ndGgodW5pcXVlKGhpZ2hfbG93X0hFSSRpZCkpKSkgKwogIHRoZW1lX21pbmltYWwoKSArCiAgdGhlbWUobGVnZW5kLnRpdGxlID0gZWxlbWVudF9ibGFuaygpKQpgYGAKCiMjIEZhdm91cml0ZSBob3VyIGFuZCBkYXkKCmBgYHtyfQpmYXZvdXJpdGVfZGF5X2hlaSA8LSBudW1iZXJfcG9zdHNfZGF5cyAlPiUKICBncm91cF9ieShpZCkgJT4lCiAgdG9wX24oMSwgY291bnQpICU+JQogIGFycmFuZ2UoaWQpCgpwcmludChmYXZvdXJpdGVfZGF5X2hlaSkKYGBgCgpgYGB7cn0KIyBDYWxjdWxhdGluZyB0aGUgbnVtYmVyIG9mIHBvc3RzIHBlciBob3VyCm51bWJlcl9wb3N0c19ob3VycyA8LSBkYXRhX3Bvc3RzX2RheXMgJT4lCiAgZ3JvdXBfYnkoaWQsIGNyZWF0ZWRfaG91cikgJT4lCiAgc3VtbWFyaXNlKGNvdW50ID0gbigpKSAlPiUKICB1bmdyb3VwKCkKCiMgSWRlbnRpZnlpbmcgdGhlIGZhdm9yaXRlIGhvdXIgZm9yIGVhY2ggSEVJCmZhdm91cml0ZV9ob3VyX2hlaSA8LSBudW1iZXJfcG9zdHNfaG91cnMgJT4lCiAgZ3JvdXBfYnkoaWQpICU+JQogIHRvcF9uKDEsIGNvdW50KSAlPiUKICBhcnJhbmdlKGlkKQoKIyBBZGRpbmcgbmV3IGNvbHVtbnMgdG8gY2F0ZWdvcml6ZSB0aGUgaG91cnMgYW5kIGFzc2lnbiBudW1lcmljYWwgdmFsdWUgZm9yIGNsdXN0ZXJpbmcKZmF2b3VyaXRlX2hvdXJfaGVpIDwtIGZhdm91cml0ZV9ob3VyX2hlaSAlPiUKICBtdXRhdGUodGltZV9vZl9kYXkgPSBjYXNlX3doZW4oCiAgICBjcmVhdGVkX2hvdXIgPj0gNSAmIGNyZWF0ZWRfaG91ciA8IDEyIH4gIk1vcm5pbmciLAogICAgY3JlYXRlZF9ob3VyID49IDEyICYgY3JlYXRlZF9ob3VyIDwgMTggfiAiQWZ0ZXJub29uIiwKICAgIFRSVUUgfiAiTmlnaHQiCiAgKSwKICB0aW1lX29mX2RheV92YWx1ZSA9IGNhc2Vfd2hlbigKICAgIGNyZWF0ZWRfaG91ciA+PSA1ICYgY3JlYXRlZF9ob3VyIDwgMTIgfiAxLAogICAgY3JlYXRlZF9ob3VyID49IDEyICYgY3JlYXRlZF9ob3VyIDwgMTggfiAzLAogICAgVFJVRSB+IDUKICApKQoKcHJpbnQoZmF2b3VyaXRlX2hvdXJfaGVpKQpgYGAKCiMjIEhlYXRtYXBzCgojIyMgRnVuY3Rpb24gdG8gcGxvdCBoZWF0bWFwIGZvciB2YXJpb3VzIEhFSXMKCmBgYHtyfQpoZWF0bWFwX21ha2VyIDwtIGZ1bmN0aW9uKHRhcmdldF9pZCl7CiAgIyBGaWx0ZXJpbmcgZGF0YSBmb3IgdGhlIHNwZWNpZmljIEhFSQogIHRhcmdldF9kYXRhIDwtIGRhdGFfcG9zdHNfZGF5cyAlPiUKICAgIGZpbHRlcihpZCA9PSB0YXJnZXRfaWQpCiAgCiAgIyBHcm91cGluZyBieSBkYXkgb2YgdGhlIHdlZWsgYW5kIGhvdXIsIGFuZCBjb3VudGluZyB0aGUgbnVtYmVyIG9mIHR3ZWV0cwogIHR3ZWV0X2NvdW50cyA8LSB0YXJnZXRfZGF0YSAlPiUKICAgIGdyb3VwX2J5KGRheV9vZl93ZWVrLCBjcmVhdGVkX2hvdXIpICU+JQogICAgc3VtbWFyaXNlKG51bV9wb3N0cyA9IG4oKSkKICAKICAjIFBsb3R0aW5nIGhlYXRtYXAKICBnZ3Bsb3QodHdlZXRfY291bnRzLCBhZXMoeCA9IGRheV9vZl93ZWVrLCB5ID0gY3JlYXRlZF9ob3VyLCBmaWxsID0gbnVtX3Bvc3RzKSkgKwogICAgZ2VvbV90aWxlKCkgKwogICAgc2NhbGVfZmlsbF9ncmFkaWVudChsb3cgPSAid2hpdGUiLCBoaWdoID0gImJsdWUiKSArCiAgICBsYWJzKHRpdGxlID0gcGFzdGUoIlBvc3QgSGVhdG1hcCBmb3IiLCB0YXJnZXRfaWQpLAogICAgICAgICB4ID0gIkRheSBvZiB0aGUgd2VlayIsCiAgICAgICAgIHkgPSAiSG91ciBvZiB0aGUgZGF5IikKfQpgYGAKCiMjIyBQbG90IG9mIGhlYXRtYXAgZm9yIGVhY2ggSEVJCgpgYGB7cn0KaGVhdG1hcF9tYWtlcigiZHVrZS5jc3YiKQpoZWF0bWFwX21ha2VyKCJlcGZsLmNzdiIpCmhlYXRtYXBfbWFrZXIoImdvZS5jc3YiKQpoZWF0bWFwX21ha2VyKCJoYXJ2YXJkLmNzdiIpCmhlYXRtYXBfbWFrZXIoImxlaWNlc3Rlci5jc3YiKQpoZWF0bWFwX21ha2VyKCJtYW5jaGVzdGVyLmNzdiIpCmhlYXRtYXBfbWFrZXIoIm1pdC5jc3YiKQpoZWF0bWFwX21ha2VyKCJzYi5jc3YiKQpoZWF0bWFwX21ha2VyKCJzdGFuZm9yZC5jc3YiKQpoZWF0bWFwX21ha2VyKCJ0cmluaXR5LmNzdiIpCmhlYXRtYXBfbWFrZXIoInd2LmNzdiIpCmhlYXRtYXBfbWFrZXIoInlhbGUuY3N2IikKYGBgCgojIyBIYXNodGFncwoKYGBge3J9CiMgVHJhbnNmb3JtaW5nIGVtcHR5IHN0cmluZ3MgaW50byBOQQpkYXRhX3Bvc3RzJGhhc2h0YWdzW2RhdGFfcG9zdHMkaGFzaHRhZ3MgPT0gIiJdIDwtIE5BCgojIFRhYmxlIHdpdGggbnVtYmVyIG9mIHVuaXF1ZSBoYXNodGFncyBhbmQgcGVyY2VudGFnZSBvZiB1c2FnZQpoYXNodGFncyA8LSBkYXRhX3Bvc3RzICU+JQogICAgICAgICAgICAgICAgZ3JvdXBfYnkoaWQpICU+JQogICAgICAgICAgICAgICAgc3VtbWFyaXNlKGNvdW50ID0gbigpLAogICAgICAgICAgICAgICAgICAgICAgICAgIG5hID0gc3VtKGlzLm5hKGhhc2h0YWdzKSksCiAgICAgICAgICAgICAgICAgICAgICAgICAgdW5pcXVlX2hhc2h0YWdzID0gbGVuZ3RoKHVuaXF1ZShoYXNodGFncykpLAogICAgICAgICAgICAgICAgICAgICAgICAgIGhhc2h0YWdfcGVyY2VudGFnZSA9IHJvdW5kKCgoY291bnQgLSBuYSkgLyBjb3VudCAqIDEwMCksIDIpKQoKcHJpbnQoaGFzaHRhZ3MpCmBgYAoKIyMjIFBsb3QgZm9yIHRoZSBjb3VudCBvZiB1bmlxdWUgaGFzaHRhZ3MgZm9yIGVhY2ggSEVJCgpgYGB7cn0KYmFycGxvdChoYXNodGFncyR1bmlxdWVfaGFzaHRhZ3MsCiAgICAgICAgbmFtZXMuYXJnID0gaGFzaHRhZ3MkaWQsCiAgICAgICAgbWFpbiA9ICJVbmlxdWUgSGFzaHRhZ3MgZm9yIEVhY2ggSEVJIiwKICAgICAgICB4bGFiID0gIkhFSSIsCiAgICAgICAgeWxhYiA9ICJDb3VudCBvZiBVbmlxdWUgSGFzaHRhZ3MiLAogICAgICAgIHlsaW0gPSBjKDAsIG1heChoYXNodGFncyR1bmlxdWVfaGFzaHRhZ3MpICsgNTApLAogICAgICAgIGxhcyA9IDIsCiAgICAgICAgY29sPSAiIzE2QTA4NSIpCgp0ZXh0KHggPSBiYXJwbG90KGhhc2h0YWdzJHVuaXF1ZV9oYXNodGFncywgcGxvdCA9IEZBTFNFKSwKICAgICB5ID0gaGFzaHRhZ3MkdW5pcXVlX2hhc2h0YWdzLAogICAgIGxhYmVscyA9IHJvdW5kKGhhc2h0YWdzJHVuaXF1ZV9oYXNodGFncywgMiksCiAgICAgcG9zID0gMykKYGBgCgojIyMgUGxvdCBmb3IgdGhlIHVzYWdlIG9mIGhhc2h0YWcgZm9yIGVhY2ggSEVJCgpgYGB7cn0KYmFycGxvdChoYXNodGFncyRoYXNodGFnX3BlcmNlbnRhZ2UsCiAgICAgICAgbmFtZXMuYXJnID0gaGFzaHRhZ3MkaWQsCiAgICAgICAgbWFpbiA9ICJIYXNodGFncyBQZXJjZW50YWdlIGZvciBFYWNoIEhFSSIsCiAgICAgICAgeGxhYiA9ICJIRUkiLAogICAgICAgIHlsYWIgPSAiSGFzaHRhZ3MgUGVyY2VudGFnZSIsCiAgICAgICAgeWxpbSA9IGMoMCwgbWF4KGhhc2h0YWdzJGhhc2h0YWdfcGVyY2VudGFnZSkgKyAzMCksCiAgICAgICAgbGFzID0gMiwKICAgICAgICBjb2w9ICIjRjFDNDBGIikKCnRleHQoeCA9IGJhcnBsb3QoaGFzaHRhZ3MkaGFzaHRhZ19wZXJjZW50YWdlLCBwbG90ID0gRkFMU0UpLAogICAgIHkgPSBoYXNodGFncyRoYXNodGFnX3BlcmNlbnRhZ2UsCiAgICAgbGFiZWxzID0gcm91bmQoaGFzaHRhZ3MkaGFzaHRhZ19wZXJjZW50YWdlLCAyKSwKICAgICBwb3MgPSAzKQpgYGAKCiMjIFVSTCB1c2FnZQoKYGBge3J9CiMgVHJhbnNmb3JtaW5nIGVtcHR5IHN0cmluZ3MgaW50byBOQQpkYXRhX3Bvc3RzJHVybHNbZGF0YV9wb3N0cyR1cmxzID09ICIiXSA8LSBOQQoKIyBUYWJsZSB3aXRoIG51bWJlciBvZiBwb3N0LCBudW1iZXIgb2YgTkEgYW5kIHVybCBwZXJjZW50YWdlIG9mIHVzYWdlCnVybF91c2FnZSA8LSBkYXRhX3Bvc3RzICU+JQogICAgICAgICAgICAgICAgZ3JvdXBfYnkoaWQpICU+JQogICAgICAgICAgICAgICAgc3VtbWFyaXNlKGNvdW50ID0gbigpLAogICAgICAgICAgICAgICAgICAgICAgICAgIG5hID0gc3VtKGlzLm5hKHVybHMpKSwKICAgICAgICAgICAgICAgICAgICAgICAgICB1cmxfcGVyY2VudGFnZSA9IHJvdW5kKCgoY291bnQgLSBuYSkgLyBjb3VudCAqIDEwMCksIDIpKQoKcHJpbnQodXJsX3VzYWdlKQpgYGAKCiMjIyBQbG90IGZvciB0aGUgdXNhZ2Ugb2YgaGFzaHRhZyBmb3IgZWFjaCBIRUkKCmBgYHtyfQpiYXJwbG90KHVybF91c2FnZSR1cmxfcGVyY2VudGFnZSwKICAgICAgICBuYW1lcy5hcmcgPSB1cmxfdXNhZ2UkaWQsCiAgICAgICAgbWFpbiA9ICJVcmxzIFBlcmNlbnRhZ2UgZm9yIEVhY2ggSEVJIiwKICAgICAgICB4bGFiID0gIkhFSSIsCiAgICAgICAgeWxhYiA9ICJVcmxzIFBlcmNlbnRhZ2UiLAogICAgICAgIHlsaW0gPSBjKDAsIG1heCh1cmxfdXNhZ2UkdXJsX3BlcmNlbnRhZ2UpICsgMTApLAogICAgICAgIGxhcyA9IDIsCiAgICAgICAgY29sPSAiIzhFNDRBRCIpCgp0ZXh0KHggPSBiYXJwbG90KHVybF91c2FnZSR1cmxfcGVyY2VudGFnZSwgcGxvdCA9IEZBTFNFKSwKICAgICB5ID0gdXJsX3VzYWdlJHVybF9wZXJjZW50YWdlLAogICAgIGxhYmVscyA9IHJvdW5kKHVybF91c2FnZSR1cmxfcGVyY2VudGFnZSwgMiksCiAgICAgcG9zID0gMykKYGBgCgojIyBUZXh0CgpgYGB7cn0KZGF0YV9wb3N0c19jb250ZW50IDwtIGRhdGFfcG9zdHMgJT4lCiAgICAgICAgICAgIHNlbGVjdChpZCwgdGV4dCkKCiMgQ291bnRpbmcgbnVtYmVyIG9mIHdvcmRzCmRhdGFfcG9zdHNfY29udGVudCA8LSBkYXRhX3Bvc3RzX2NvbnRlbnQgJT4lCiAgbXV0YXRlKG51bV93b3JkcyA9IGxlbmd0aHMoc3Ryc3BsaXQodGV4dCwgIlxccysiKSkpCgojIEdyb3VwaW5nIGJ5IEhFSSBhbmQgY2FsY3VsYXRlIGF2ZXJhZ2UsIG1pbmltdW0sIGFuZCBtYXhpbXVtIHZhbHVlcyBvZiBudW1iZXIgb2Ygd29yZHMKZGF0YV9wb3N0c19jb250ZW50X21ldHJpY3MgPC0gZGF0YV9wb3N0c19jb250ZW50ICU+JQogIGdyb3VwX2J5KGlkKSAlPiUKICBzdW1tYXJpc2UoYXZlcmFnZV9udW1fd29yZHMgPSBtZWFuKG51bV93b3JkcyksCiAgICAgICAgICAgIG1pbl9udW1fd29yZHMgPSBtaW4obnVtX3dvcmRzKSwKICAgICAgICAgICAgbWF4X251bV93b3JkcyA9IG1heChudW1fd29yZHMpKQpwcmludChkYXRhX3Bvc3RzX2NvbnRlbnRfbWV0cmljcykKYGBgCgojIyMgUGxvdCBmb3IgdGhlIGF2ZXJhZ2UsIG1heGltdW0gYW5kIG1pbmltdW0gdmFsdWVzIG9mIHdvcmRzIGZvciBlYWNoIEhFSQoKYGBge3J9CmdncGxvdChkYXRhX3Bvc3RzX2NvbnRlbnRfbWV0cmljcywgYWVzKHggPSBpZCwgeSA9IGF2ZXJhZ2VfbnVtX3dvcmRzKSkgKwogIGdlb21fcG9pbnQoYWVzKGNvbG9yID0gIkF2ZXJhZ2UiKSkgKwogIGdlb21fZXJyb3JiYXIoYWVzKHltaW4gPSBtaW5fbnVtX3dvcmRzLCB5bWF4ID0gbWF4X251bV93b3JkcywgY29sb3IgPSAiUmFuZ2UiKSwgd2lkdGggPSAwLjIpICsKICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gYygiQXZlcmFnZSIgPSAiIzE5NzZEMiIsICJSYW5nZSIgPSAiI0VGNTM1MCIpKSArCiAgbGFicyh0aXRsZSA9ICJXb3JkIENvdW50IFN1bW1hcnkgYnkgSEVJIiwKICAgICAgIHggPSAiSEVJIiwKICAgICAgIHkgPSAiTnVtYmVyIG9mIFdvcmRzIiwKICAgICAgIGNvbG9yID0gIk1ldHJpYyIpICsKICB0aGVtZV9taW5pbWFsKCkgKwogIHRoZW1lKGF4aXMudGV4dC54ID0gZWxlbWVudF90ZXh0KGFuZ2xlID0gOTAsIHZqdXN0ID0gMC41LCBoanVzdCA9IDEpKQpgYGAKCiMgUG9zdHMgQ2xhc3NpZmljYXRpb24KCiMjIEZ1bmN0aW9uIHRvIGNsZWFudXAgdGV4dAoKYGBge3J9CmNsZWFudXAgPC0gZnVuY3Rpb24oZG9jcywgc3BlYy53b3Jkcz1OVUxMKXsKICAjIGxvd2VyY2FzZQogIGRvY3MgPC0gdG1fbWFwKGRvY3MsIGNvbnRlbnRfdHJhbnNmb3JtZXIodG9sb3dlcikpCiAgIyBybSBudW1iZXJzCiAgZG9jcyA8LSB0bV9tYXAoZG9jcywgcmVtb3ZlTnVtYmVycykKICAjIHJtIGVuZ2xpc2ggY29tbW9uIHN0b3BXb3JkcwogIGRvY3MgPC0gdG1fbWFwKGRvY3MsIHJlbW92ZVdvcmRzLCBzdG9wd29yZHMoImVuZ2xpc2giKSkKICAjIGlmIHN0b3B3b3JkcyBhcmUgc3BlY2lmaWVkIGFzIGEgY2hhcmFjdGVyIHZlY3RvcgogIGlmKCFpcy5udWxsKHNwZWMud29yZHMpKQogICAgZG9jcyA8LSB0bV9tYXAoZG9jcywgcmVtb3ZlV29yZHMsIHNwZWMud29yZHMpCiAgIyBybSBwdW5jdHVhdGlvbnMKICBkb2NzIDwtIHRtX21hcChkb2NzLCByZW1vdmVQdW5jdHVhdGlvbikKICAjIHJtIGV4dHJhIHdoaXRlIHNwYWNlcwogIGRvY3MgPC0gdG1fbWFwKGRvY3MsIHN0cmlwV2hpdGVzcGFjZSkKICAjIGxlbW1hdGl6aW5nIHRleHQKICBkb2NzIDwtIHRtX21hcChkb2NzLCBsZW1tYXRpemVfd29yZHMpCiAgCiAgZG9jcwp9CmBgYAoKIyMgRnVjbnRpb24gdG8gcGFpciB3b3JkcwoKYGBge3J9CnBhaXJpbmcgPC0gZnVuY3Rpb24od29yZF9kaWN0LCBwYWlyZWRfd29yZHMpIHsKICBpZiAobnJvdyh3b3JkX2RpY3QpICE9IGxlbmd0aChwYWlyZWRfd29yZHMpKSB7CiAgICBzdG9wKCJUaGUgbnVtYmVyIG9mIHJvd3MgaW4gd29yZF9kaWN0IGFuZCB0aGUgbGVuZ3RoIG9mIHBhaXJlZF93b3JkcyBzaG91bGQgYmUgZXF1YWwuIikKICB9CiAgCiAgd29yZF9kaWN0IDwtIHRpYmJsZTo6dGliYmxlKHdvcmQgPSB3b3JkX2RpY3QkdmFsdWUpCiAgcmVzdWx0IDwtIHRpYmJsZTo6dGliYmxlKHdvcmQgPSB3b3JkX2RpY3Qkd29yZCwgY2F0ZWdvcnkgPSBwYWlyZWRfd29yZHMpCiAgcmV0dXJuKHJlc3VsdCkKfQpgYGAKCiMjIEZ1bmN0aW9uIHRvIGNsYXNzaWZ5IHRleHQKCmBgYHtyfQpjbGFzc2lmeV90ZXh0IDwtIGZ1bmN0aW9uKHRleHQsIHdvcmRfcGFpcikgewogICMgVG9rZW5pemluZyB0aGUgdGV4dCBhbmQgY29udmVydGluZyB0byBsb3dlcmNhc2UKICB3b3JkcyA8LSB0b2xvd2VyKHVubGlzdChzdHJzcGxpdCh0ZXh0LCAiXFxXKyIpKSkKICAKICAjIEZpbmRpbmcgdGhlIGZyZXF1ZW50IHdvcmRzIGluIHRoZSB0ZXh0CiAgZnJlcV93b3JkcyA8LSB3b3Jkc1t3b3JkcyAlaW4lIHdvcmRfcGFpciR3b3JkXQogIAogIGlmIChsZW5ndGgoZnJlcV93b3JkcykgPT0gMCkgewogICAgcmV0dXJuKCJVbmtub3duIikgICMgUmV0dXJuaW5nIFVua25vd24gaWYgbm8gZnJlcXVlbnQgd29yZHMgYXJlIGZvdW5kCiAgfQogIAogICMgR2V0dGluZyB0aGUgY29ycmVzcG9uZGluZyBjYXRlZ29yaWVzCiAgY2F0ZWdvcmllcyA8LSB3b3JkX3BhaXIkY2F0ZWdvcnlbd29yZF9wYWlyJHdvcmQgJWluJSBmcmVxX3dvcmRzXQogIAogICMgU29ydGluZyBjYXRlZ29yaWVzIGJ5IGZyZXF1ZW5jeSBpbiBkZXNjZW5kaW5nIG9yZGVyIGFuZCByZXR1cm4gdGhlIG1vc3QgZnJlcXVlbnQgb25lCiAgcmV0dXJuKG5hbWVzKHNvcnQodGFibGUoY2F0ZWdvcmllcyksIGRlY3JlYXNpbmcgPSBUUlVFKSlbMV0pCn0KYGBgCgojIyBGdW5jdGlvbiB0byBjcmVhdGUgY29ycHVzCgpgYGB7cn0KY29ycHVzX21ha2VyIDwtIGZ1bmN0aW9uKHRleHQpIHsKICB0ZXh0cyA8LSB0ZXh0JHRleHQKICB2YyA8LSBWZWN0b3JTb3VyY2UodGV4dHMpCiAgY29ycHVzIDwtIENvcnB1cyh2YykKICAKICByZXR1cm4oY29ycHVzKQp9CmBgYAoKIyMgRnVuY3Rpb24gdG8gZmluZCBmcmVxdWVudCB0ZXJtcyBpbiB0ZXh0cwoKYGBge3J9CmZyZXFfdGVybXMgPC0gZnVuY3Rpb24oY2xlYW5fdGV4dCwgbnVtYmVyKSB7CiAgZHRtIDwtIERvY3VtZW50VGVybU1hdHJpeChjbGVhbl90ZXh0KQogIGR0bS50ZmlkZiA8LSB3ZWlnaHRUZklkZihkdG0pCiAgCiAgbWRmIDwtIGFzX3RpYmJsZShhcy5tYXRyaXgoZHRtLnRmaWRmKSkKICAKICBtZGYuZnJlcSA8LSBtZGYgJT4lCiAgc2VsZWN0KGZpbmRGcmVxVGVybXMoZHRtLCBudW1iZXIpKSAlPiUKICBzdW1tYXJpc2VfYWxsKHN1bSkgJT4lCiAgZ2F0aGVyKCkgJT4lCiAgYXJyYW5nZShkZXNjKHZhbHVlKSkKCiAgbWRmLmZyZXEka2V5IDwtIAogICAgZmFjdG9yKG1kZi5mcmVxJGtleSwKICAgICAgICAgICBsZXZlbHMgPSBtZGYuZnJlcSRrZXlbb3JkZXIobWRmLmZyZXEkdmFsdWUpXSkKICAKICB3b3JkX2RpY3Rpb25hcnkgPC0gYXNfdGliYmxlKG1kZi5mcmVxJGtleSkKICAKICBnZ3Bsb3QobWRmLmZyZXEsIGFlcyh4PWtleSwgeT12YWx1ZSkpICsKICAgIGdlb21fYmFyKHN0YXQ9ImlkZW50aXR5IikgKwogICAgbGFicyh4PSJ0ZXJtcyIsIHk9ImZyZXEiKSArIGNvb3JkX2ZsaXAoKQogIAogIHJldHVybih3b3JkX2RpY3Rpb25hcnkpCn0KYGBgCgojIyMgRHVrZQoKYGBge3J9CmR1a2VfdGV4dCA8LSBzdWJzZXQoZGF0YV9wb3N0c19jb250ZW50LCBpZCA9PSAiZHVrZS5jc3YiKQpkdWtlX2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIoZHVrZV90ZXh0KQoKZHVrZV9jbGVhbl90ZXh0IDwtY2xlYW51cChkdWtlX2NvcnB1cywgYygibmV3IiwgIndpbGwiLCAiY2hhbmdlIiwgIm5vcnRoIiwgImNhbiIsICJmaXJzdCIsICJ5ZWFyIiwgImNhcm9saW5hIiwgInllYXJzIiwgInN1bW1lciIsICJtZWV0IiwgIm9uZSIpKQoKZHVrZV93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhkdWtlX2NsZWFuX3RleHQsIDMwKQoKZHVrZV93b3JkX2NhdGVnb3J5IDwtIGMoIkltYWdlIiwgIkVkdWNhdGlvbiIsICJFZHVjYXRpb24iLCAiSW1hZ2UiLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCAiUmVzZWFyY2giLCAiSW1hZ2UiLCAiSW1hZ2UiLCAiRW5nYWdlbWVudCIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJJbWFnZSIsICJFZHVjYXRpb24iLCAiRWR1Y2F0aW9uIiwgIlNvY2lldHkiLCAiRWR1Y2F0aW9uIiwgIkVuZ2FnZW1lbnQiKQoKZHVrZV93b3JkX3BhaXIgPC0gcGFpcmluZyhkdWtlX3dvcmRfZGljdGlvbmFyeSwgZHVrZV93b3JkX2NhdGVnb3J5KQoKcHJpbnQoZHVrZV93b3JkX3BhaXIpCgojIEFwcGx5aW5nIHRoZSBjbGFzc2lmeV90ZXh0IGZ1bmN0aW9uIHRvIGVhY2ggdGV4dApkdWtlX3RleHQgPC0gZHVrZV90ZXh0ICU+JQogIG11dGF0ZShjYXRlZ29yeSA9IHNhcHBseSh0ZXh0LCBjbGFzc2lmeV90ZXh0LCB3b3JkX3BhaXIgPSBkdWtlX3dvcmRfcGFpcikpCgpwcmludChkdWtlX3RleHQpCmBgYAoKIyMjIGVwZmwKCmBgYHtyfQplcGZsX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gImVwZmwuY3N2IikKZXBmbF9jb3JwdXMgPC0gY29ycHVzX21ha2VyKGVwZmxfdGV4dCkKCmVwZmxfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKGVwZmxfY29ycHVzLCBjKCJuZXciLCAiYW1wIiwgImNhbiIsICJ3aWxsIiwgIm9uZSIsICJub3ciLCAi4oCTIiwgInZpYSIsICJwb3J0cmFpdCIpKQoKZXBmbF93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhlcGZsX2NsZWFuX3RleHQsIDMwKQoKZXBmbF93b3JkX2NhdGVnb3J5IDwtIGMoIkltYWdlIiwgIlJlc2VhcmNoIiwgIkltYWdlIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJJbWFnZSIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iLCAiUmVzZWFyY2giLCBOQSwgIkVkdWNhdGlvbiIsICJFZHVjYXRpb24iKQoKZXBmbF93b3JkX3BhaXIgPC0gcGFpcmluZyhlcGZsX3dvcmRfZGljdGlvbmFyeSwgZXBmbF93b3JkX2NhdGVnb3J5KQoKcHJpbnQoZXBmbF93b3JkX3BhaXIpCgojIEFwcGx5aW5nIHRoZSBjbGFzc2lmeV90ZXh0IGZ1bmN0aW9uIHRvIGVhY2ggdGV4dAplcGZsX3RleHQgPC0gZXBmbF90ZXh0ICU+JQogIG11dGF0ZShjYXRlZ29yeSA9IHNhcHBseSh0ZXh0LCBjbGFzc2lmeV90ZXh0LCB3b3JkX3BhaXIgPSBlcGZsX3dvcmRfcGFpcikpCgpwcmludChlcGZsX3RleHQpCmBgYAoKIyMjIGdvZQoKYGBge3J9CmdvZV90ZXh0IDwtIHN1YnNldChkYXRhX3Bvc3RzX2NvbnRlbnQsIGlkID09ICJnb2UuY3N2IikKZ29lX2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIoZ29lX3RleHQpCgpnb2VfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKGdvZV9jb3JwdXMsIGMoImNhbiIsICJuZXciLCAiYW1wIiwgIndpbGwiKSkKCmdvZV93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhnb2VfY2xlYW5fdGV4dCwgMTUpCgpnb2Vfd29yZF9jYXRlZ29yeSA8LSBjKCJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJJbWFnZSIsICJFZHVjYXRpb24iLCAiRW5nYWdlbWVudCIsICJJbWFnZSIsICJFbmdhZ2VtZW50IiwgIkVkdWNhdGlvbiIpCgpnb2Vfd29yZF9wYWlyIDwtIHBhaXJpbmcoZ29lX3dvcmRfZGljdGlvbmFyeSwgZ29lX3dvcmRfY2F0ZWdvcnkpCgpwcmludChnb2Vfd29yZF9wYWlyKQoKIyBBcHBseWluZyB0aGUgY2xhc3NpZnlfdGV4dCBmdW5jdGlvbiB0byBlYWNoIHRleHQKZ29lX3RleHQgPC0gZ29lX3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IGdvZV93b3JkX3BhaXIpKQoKcHJpbnQoZ29lX3RleHQpCmBgYAoKIyMjIGhhcnZhcmQKCmBgYHtyfQpoYXJ2YXJkX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gImhhcnZhcmQuY3N2IikKaGFydmFyZF9jb3JwdXMgPC0gY29ycHVzX21ha2VyKGhhcnZhcmRfdGV4dCkKCmhhcnZhcmRfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKGhhcnZhcmRfY29ycHVzLCBjKCJuZXciLCAiY2FuIiwgIndpbGwiLCAic3VtbWVyIiwgInllYXIiLCAiZmlyc3QiLCAibWF5IiwgIi0iLCAicmVjZW50IiwgInllYXJzIiwgIm9uZSIsICJzYWlkIiwgInRpbWUiLCAibWFueSIsICJ3b3JsZCIsICJjaGFuZ2UiKSkKCmhhcnZhcmRfd29yZF9kaWN0aW9uYXJ5IDwtIGZyZXFfdGVybXMoaGFydmFyZF9jbGVhbl90ZXh0LCA1MCkKCmhhcnZhcmRfd29yZF9jYXRlZ29yeSA8LSBjKCJJbWFnZSIsICJFZHVjYXRpb24iLCAiUmVzZWFyY2giLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCAiUmVzZWFyY2giLCAiRW5nYWdlbWVudCIsICJFZHVjYXRpb24iLCAiU29jaWV0eSIsICJFZHVjYXRpb24iLCAiU29jaWV0eSIsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCBOQSwgIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJFZHVjYXRpb24iLCBOQSwgIlJlc2VhcmNoIiwgIlNvY2lldHkiLCAiU29jaWV0eSIsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iKQoKaGFydmFyZF93b3JkX3BhaXIgPC0gcGFpcmluZyhoYXJ2YXJkX3dvcmRfZGljdGlvbmFyeSwgaGFydmFyZF93b3JkX2NhdGVnb3J5KQoKcHJpbnQoaGFydmFyZF93b3JkX3BhaXIpCgpoYXJ2YXJkX3RleHQgPC0gaGFydmFyZF90ZXh0ICU+JQogIG11dGF0ZShjYXRlZ29yeSA9IHNhcHBseSh0ZXh0LCBjbGFzc2lmeV90ZXh0LCB3b3JkX3BhaXIgPSBoYXJ2YXJkX3dvcmRfcGFpcikpCgpwcmludChoYXJ2YXJkX3RleHQpCmBgYAoKIyMjIGxlaWNlc3RlcgoKYGBge3J9CmxlaWNlc3Rlcl90ZXh0IDwtIHN1YnNldChkYXRhX3Bvc3RzX2NvbnRlbnQsIGlkID09ICJsZWljZXN0ZXIuY3N2IikKbGVpY2VzdGVyX2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIobGVpY2VzdGVyX3RleHQpCgpsZWljZXN0ZXJfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKGxlaWNlc3Rlcl9jb3JwdXMsIGMoIvCfkYkiLCAiZGF5IiwgIm5ldyIsICJjbGVhciIsICJ5ZWFyIiwgImNhbiIsICJ3aWxsIiwgInRpbWUiLCAic3BhY2UiLCAib25lIiwgImZpcnN0IikpCgpsZWljZXN0ZXJfd29yZF9kaWN0aW9uYXJ5IDwtIGZyZXFfdGVybXMobGVpY2VzdGVyX2NsZWFuX3RleHQsIDYwKQoKbGVpY2VzdGVyX3dvcmRfY2F0ZWdvcnkgPC0gYygiSW1hZ2UiLCBOQSwgIkltYWdlIiwgIkVuZ2FnZW1lbnQiLCAiU29jaWV0eSIsICJFZHVjYXRpb24iLCAiRWR1Y2F0aW9uIiwgIlJlc2VhcmNoIiwgIkVuZ2FnZW1lbnQiLCAiRW5nYWdlbWVudCIsICJJbWFnZSIsICJFbmdhZ2VtZW50IiwgIkVkdWNhdGlvbiIsICJFbmdhbWVtZW50IiwgIkVuZ2FnZW1lbnQiLCAiRWR1Y2F0aW9uIiwgIkltYWdlIiwgIkVuZ2FnZW1lbnQiLCAiRWR1Y2F0aW9uIiwgTkEsICJJbWFnZSIpCgpsZWljZXN0ZXJfd29yZF9wYWlyIDwtIHBhaXJpbmcobGVpY2VzdGVyX3dvcmRfZGljdGlvbmFyeSwgbGVpY2VzdGVyX3dvcmRfY2F0ZWdvcnkpCgpwcmludChsZWljZXN0ZXJfd29yZF9wYWlyKQoKbGVpY2VzdGVyX3RleHQgPC0gbGVpY2VzdGVyX3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IGxlaWNlc3Rlcl93b3JkX3BhaXIpKQoKcHJpbnQobGVpY2VzdGVyX3RleHQpCmBgYAoKIyMjIG1hbmNoZXN0ZXIKCmBgYHtyfQptYW5jaGVzdGVyX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gIm1hbmNoZXN0ZXIuY3N2IikKbWFuY2hlc3Rlcl9jb3JwdXMgPC0gY29ycHVzX21ha2VyKG1hbmNoZXN0ZXJfdGV4dCkKCm1hbmNoZXN0ZXJfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKG1hbmNoZXN0ZXJfY29ycHVzLCBjKCLwn5GHIiwgImNhbiIsICJqdXN0IiwgIndpbGwiLCAiZ2V0IiwgIndlbGwiLCAib25lIiwgImhlbHAiLCAibm93IiwgIm5ldyIsICJyZWFkIiwgImNvbmdyYXR1bGF0aW9ucyIpKQoKbWFuY2hlc3Rlcl93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhtYW5jaGVzdGVyX2NsZWFuX3RleHQsIDUwKQoKbWFuY2hlc3Rlcl93b3JkX2NhdGVnb3J5IDwtIGMoTkEsICJJbWFnZSIsICJFbmdhZ2VtZW50IiwgIkltYWdlIiwgIkVkdWNhdGlvbiIsICJFZHVjYXRpb24iLCAiRW5nYWdlbWVudCIsICJJbWFnZSIsICJFZHVjYXRpb24iLCAiU29jaWV0eSIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iLCAiRWR1Y2F0aW9uIiwgIkVkdWNhdGlvbiIsICJFbmdhZ2VtZW50IiwgIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJSZXNlYXJjaCIpCgptYW5jaGVzdGVyX3dvcmRfcGFpciA8LSBwYWlyaW5nKG1hbmNoZXN0ZXJfd29yZF9kaWN0aW9uYXJ5LCBtYW5jaGVzdGVyX3dvcmRfY2F0ZWdvcnkpCgpwcmludChtYW5jaGVzdGVyX3dvcmRfcGFpcikKCm1hbmNoZXN0ZXJfdGV4dCA8LSBtYW5jaGVzdGVyX3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IG1hbmNoZXN0ZXJfd29yZF9wYWlyKSkKCnByaW50KG1hbmNoZXN0ZXJfdGV4dCkKYGBgCgojIyMgbWl0CgpgYGB7cn0KbWl0X3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gIm1pdC5jc3YiKQptaXRfY29ycHVzIDwtIGNvcnB1c19tYWtlcihtaXRfdGV4dCkKCm1pdF9jbGVhbl90ZXh0IDwtIGNsZWFudXAobWl0X2NvcnB1cywgYygibmV3IiwgImNhbiIsICJzYXlzIiwgIuKAnCIsICLigJkiLCAi4oCUIiwgIiIsICJtYXkiLCAiZmlyc3QiLCAid2lsbCIsICJ1c2luZyIsICJ3YXkiLCAib25lIiwgInNjaWVuY2UiKSkKCm1pdF93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhtaXRfY2xlYW5fdGV4dCwgNDApCgptaXRfd29yZF9jYXRlZ29yeSA8LSBjKCJJbWFnZSIsICJSZXNlYXJjaCIsIE5BLCAiSW1hZ2UiLCBOQSwgIkVkdWNhdGlvbiIsICJSZXNlYXJjaCIsIE5BLCAiRWR1Y2F0aW9uIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJSZXNlYXJjaCIsICJSZXNlYXJjaCIsICJJbWFnZSIsICJFZHVjYXRpb24iLCAiRWR1Y2F0aW9uIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlNvY2lldHkiKQoKbWl0X3dvcmRfcGFpciA8LSBwYWlyaW5nKG1pdF93b3JkX2RpY3Rpb25hcnksIG1pdF93b3JkX2NhdGVnb3J5KQoKcHJpbnQobWl0X3dvcmRfcGFpcikKCm1pdF90ZXh0IDwtIG1pdF90ZXh0ICU+JQogIG11dGF0ZShjYXRlZ29yeSA9IHNhcHBseSh0ZXh0LCBjbGFzc2lmeV90ZXh0LCB3b3JkX3BhaXIgPSBtaXRfd29yZF9wYWlyKSkKCnByaW50KG1pdF90ZXh0KQpgYGAKCiMjIyBzYgoKYGBge3J9CnNiX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gInNiLmNzdiIpCnNiX2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIoc2JfdGV4dCkKCnNiX2NsZWFuX3RleHQgPC0gY2xlYW51cChzYl9jb3JwdXMsIGMoIm5ldyIsICLigJQiLCAid2lsbCIsICJ3ZWVrIiwgImNhbiIsICJhbXAiLCAidmlhIiwgIm5vdyIsICJmdXR1cmUiLCAiZmlyc3QiKSkKCnNiX3dvcmRfZGljdGlvbmFyeSA8LSBmcmVxX3Rlcm1zKHNiX2NsZWFuX3RleHQsIDMwKQoKc2Jfd29yZF9jYXRlZ29yeSA8LSBjKCJJbWFnZSIsICJSZXNlYXJjaCIsICJJbWFnZSIsICJJbWFnZSIsICJJbWFnZSIsIE5BLCAiUmVzZWFyY2giLCAiSW1hZ2UiLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCAiRWR1Y2F0aW9uIiwgIkltYWdlIiwgIkVkdWNhdGlvbiIsICJJbWFnZSIsICJFZHVjYXRpb24iLCAiSW1hZ2UiLCAiU29jaWV0eSIpCgpzYl93b3JkX3BhaXIgPC0gcGFpcmluZyhzYl93b3JkX2RpY3Rpb25hcnksIHNiX3dvcmRfY2F0ZWdvcnkpCgpwcmludChzYl93b3JkX3BhaXIpCgpzYl90ZXh0IDwtIHNiX3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IHNiX3dvcmRfcGFpcikpCgpwcmludChzYl90ZXh0KQpgYGAKCiMjIyBzdGFuZm9yZAoKYGBge3J9CnN0YW5mb3JkX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gInN0YW5mb3JkLmNzdiIpCnN0YW5mb3JkX2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIoc3RhbmZvcmRfdGV4dCkKCnN0YW5mb3JkX2NsZWFuX3RleHQgPC0gY2xlYW51cChzdGFuZm9yZF9jb3JwdXMsIGMoIm5ldyIsICJ3aWxsIikpCgpzdGFuZm9yZF93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyhzdGFuZm9yZF9jbGVhbl90ZXh0LCAyNSkKCnN0YW5mb3JkX3dvcmRfY2F0ZWdvcnkgPC0gYygiSW1hZ2UiLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCAiRWR1Y2F0aW9uIiwgTkEpCgpzdGFuZm9yZF93b3JkX3BhaXIgPC0gcGFpcmluZyhzdGFuZm9yZF93b3JkX2RpY3Rpb25hcnksIHN0YW5mb3JkX3dvcmRfY2F0ZWdvcnkpCgpwcmludChzdGFuZm9yZF93b3JkX3BhaXIpCgpzdGFuZm9yZF90ZXh0IDwtIHN0YW5mb3JkX3RleHQgJT4lCiAgbXV0YXRlKGNhdGVnb3J5ID0gc2FwcGx5KHRleHQsIGNsYXNzaWZ5X3RleHQsIHdvcmRfcGFpciA9IHN0YW5mb3JkX3dvcmRfcGFpcikpCgpwcmludChzdGFuZm9yZF90ZXh0KQpgYGAKCiMjIyB0cmluaXR5CgpgYGB7cn0KdHJpbml0eV90ZXh0IDwtIHN1YnNldChkYXRhX3Bvc3RzX2NvbnRlbnQsIGlkID09ICJ0cmluaXR5LmNzdiIpCnRyaW5pdHlfY29ycHVzIDwtIGNvcnB1c19tYWtlcih0cmluaXR5X3RleHQpCgp0cmluaXR5X2NsZWFuX3RleHQgPC0gY2xlYW51cCh0cmluaXR5X2NvcnB1cywgYygiYW1wIiwgInJlYWQiLCAibmV3IiwgImNhbiIsICJ3aWxsIiwgIndlZWsiLCAid29yayIsICJncmVhdCIsICJkYXkiLCAidmlzaXQiLCAiaXJpc2giLCAiZmlyc3QiLCAibGVkIiwgImNvbmdyYXR1bGF0aW9ucyIpKQoKdHJpbml0eV93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyh0cmluaXR5X2NsZWFuX3RleHQsIDQwKQoKdHJpbml0eV93b3JkX2NhdGVnb3J5IDwtIGMoIkltYWdlIiwgIkVkdWNhdGlvbiIsICJSZXNlYXJjaCIsICJJbWFnZSIsICJSZXNlYXJjaCIsICJJbWFnZSIsICJSZXNlYXJjaCIsICJTb2NpZXR5IiwgIlNvY2lldHkiLCAiRWR1Y2F0aW9uIiwgIkVuZ2FnZW1lbnQiLCAiRW5nYWdlbWVudCIsICJFZHVjYXRpb24iLCAiSW1hZ2UiLCAiRWR1Y2F0aW9uIiwgIkltYWdlIiwgIlJlc2VhcmNoIikKCnRyaW5pdHlfd29yZF9wYWlyIDwtIHBhaXJpbmcodHJpbml0eV93b3JkX2RpY3Rpb25hcnksIHRyaW5pdHlfd29yZF9jYXRlZ29yeSkKCnByaW50KHRyaW5pdHlfd29yZF9wYWlyKQoKdHJpbml0eV90ZXh0IDwtIHRyaW5pdHlfdGV4dCAlPiUKICBtdXRhdGUoY2F0ZWdvcnkgPSBzYXBwbHkodGV4dCwgY2xhc3NpZnlfdGV4dCwgd29yZF9wYWlyID0gdHJpbml0eV93b3JkX3BhaXIpKQoKcHJpbnQodHJpbml0eV90ZXh0KQpgYGAKCiMjIyB3dgoKYGBge3J9Cnd2X3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gInd2LmNzdiIpCnd2X2NvcnB1cyA8LSBjb3JwdXNfbWFrZXIod3ZfdGV4dCkKCnd2X2NsZWFuX3RleHQgPC0gY2xlYW51cCh3dl9jb3JwdXMsIGMoIvCfkpvwn5KZIiwgIvCfmYwiLCAi8J+RiSIsICJoYXBweSIsICJkYXkiLCAic2VlIiwgIndlZWsiLCAiZ3JlYXQiLCAiY2FuIiwgIndpbGwiLCAid2VsbCIsICJrbm93IiwgIm5ldyIsICJub3ciLCAiZ2V0IiwgImp1c3QiKSkKCnd2X3dvcmRfZGljdGlvbmFyeSA8LSBmcmVxX3Rlcm1zKHd2X2NsZWFuX3RleHQsIDQwKQoKd3Zfd29yZF9jYXRlZ29yeSA8LSBjKE5BLCAiSW1hZ2UiLCAiSW1hZ2UiLCBOQSwgIkltYWdlIiwgIkltYWdlIiwgIkVuZ2FnZW1lbnQiLCAiRWR1Y2F0aW9uIiwgIkVkdWNhdGlvbiIsICJFZHVjYXRpb24iLCBOQSkKCnd2X3dvcmRfcGFpciA8LSBwYWlyaW5nKHd2X3dvcmRfZGljdGlvbmFyeSwgd3Zfd29yZF9jYXRlZ29yeSkKCnByaW50KHd2X3dvcmRfcGFpcikKCnd2X3RleHQgPC0gd3ZfdGV4dCAlPiUKICBtdXRhdGUoY2F0ZWdvcnkgPSBzYXBwbHkodGV4dCwgY2xhc3NpZnlfdGV4dCwgd29yZF9wYWlyID0gd3Zfd29yZF9wYWlyKSkKCnByaW50KHd2X3RleHQpCmBgYAoKIyMjIHlhbGUKCmBgYHtyfQp5YWxlX3RleHQgPC0gc3Vic2V0KGRhdGFfcG9zdHNfY29udGVudCwgaWQgPT0gInlhbGUuY3N2IikKeWFsZV9jb3JwdXMgPC0gY29ycHVzX21ha2VyKHlhbGVfdGV4dCkKCnlhbGVfY2xlYW5fdGV4dCA8LSBjbGVhbnVwKHlhbGVfY29ycHVzLCBjKCJuZXciLCAi4oCUIiwgIndpbGwiLCAiY2FuIiwgIiciLCAiZmlyc3QiLCAid29yayIsICJyZWFkIiwgImhlbHAiLCAieWVhciIpKQoKeWFsZV93b3JkX2RpY3Rpb25hcnkgPC0gZnJlcV90ZXJtcyh5YWxlX2NsZWFuX3RleHQsIDcwKQoKeWFsZV93b3JkX2NhdGVnb3J5IDwtIGMoIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJFZHVjYXRpb24iLCBOQSwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIlJlc2VhcmNoIiwgIkVkdWNhdGlvbiIsICJJbWFnZSIsICJSZXNlYXJjaCIsICJFZHVjYXRpb24iLCAiSW1hZ2UiLCAiUmVzZWFyY2giLCBOQSwgIlJlc2VhcmNoIiwgIkltYWdlIiwgIkltYWdlIiwgIlJlc2VhcmNoIiwgIkltYWdlIiwgIlJlc2VhcmNoIiwgIlNvY2lldHkiLCAiUmVzZWFyY2giLCAiU29jaWV0eSIsICJTb2NpZXR5IikKCnlhbGVfd29yZF9wYWlyIDwtIHBhaXJpbmcoeWFsZV93b3JkX2RpY3Rpb25hcnksIHlhbGVfd29yZF9jYXRlZ29yeSkKCnByaW50KHlhbGVfd29yZF9wYWlyKQoKeWFsZV90ZXh0IDwtIHlhbGVfdGV4dCAlPiUKICBtdXRhdGUoY2F0ZWdvcnkgPSBzYXBwbHkodGV4dCwgY2xhc3NpZnlfdGV4dCwgd29yZF9wYWlyID0geWFsZV93b3JkX3BhaXIpKQoKcHJpbnQoeWFsZV90ZXh0KQpgYGAKCiMgU2VudGltZW50IEFuYWx5c2lzCgojIyMgRnVuY3Rpb24gZm9yIGVtb3Rpb25zCgpgYGB7cn0KZW1vdGlvbnNfbWFrZXIgPC0gZnVuY3Rpb24odGV4dHMsIGhlaV9uYW1lKXsKICB0ZXh0X2Vtb3Rpb24gPC0gZ2V0X25yY19zZW50aW1lbnQodGV4dHMpCiAgCiAgIyBQcm9wb3J0aW9ucyBmb3IgdGV4dCB2YWx1ZXMgCiAgZW1vdGlvbl9wcm9wb3J0aW9ucyA8LSBjb2xTdW1zKHByb3AudGFibGUodGV4dF9lbW90aW9uWywgMTo4XSkpCiAgCiAgIyBWaXN1YWxpemF0aW9uIG9mIHBlcmNlbnRhZ2VzIG9uIGVhY2ggZW1vdGlvbiBmb3VuZCBvbiBwb3N0cwogIGJhcnBsb3QgPC0gYmFycGxvdCgKICAgIHNvcnQoY29sU3Vtcyhwcm9wLnRhYmxlKHRleHRfZW1vdGlvblssIDE6OF0pKSksIAogICAgaG9yaXogPSBUUlVFLCAKICAgIGNleC5uYW1lcyA9IDAuNywgCiAgICBsYXMgPSAxLCAKICAgIG1haW4gPSBzcHJpbnRmKCJFbW90aW9ucyBmb3VuZCBpbiAlcydzIHRleHRzIiwgaGVpX25hbWUpLCAKICAgIHhsYWI9IlBlcmNlbnRhZ2UiLAogICAgeGxpbSA9IGMoMCwgbWF4KGVtb3Rpb25fcHJvcG9ydGlvbnMpICogMS4yKQogICkKICAKICB0ZXh0KAogICAgeCA9IHNvcnQoZW1vdGlvbl9wcm9wb3J0aW9ucyksCiAgICB5ID0gYmFycGxvdCwgICAgICAgICAgICAgICAgICAgIAogICAgbGFiZWxzID0gc3ByaW50ZigiJS4yZiUlIiwgMTAwICogc29ydChlbW90aW9uX3Byb3BvcnRpb25zKSksCiAgICBwb3MgPSA0LCAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAKICAgIGNleCA9IDAuNyAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIAogICkKfQpgYGAKCiMjIyBEdWtlCgpgYGB7cn0KZHVrZV90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KGR1a2VfdGV4dCR0ZXh0KSwgMikKCnByaW50KGR1a2VfdGV4dCkKYGBgCgpgYGB7cn0KZW1vdGlvbnNfbWFrZXIoZHVrZV90ZXh0JHRleHQsICJkdWtlIikKYGBgCgojIyMgZXBmbAoKYGBge3J9CmVwZmxfdGV4dCRzZW50aW1lbnQgPC0gcm91bmQoZ2V0X3NlbnRpbWVudChlcGZsX3RleHQkdGV4dCksIDIpCgpwcmludChlcGZsX3RleHQpCmBgYAoKYGBge3J9CmVtb3Rpb25zX21ha2VyKGVwZmxfdGV4dCR0ZXh0LCAiZXBmbCIpCmBgYAoKIyMjIGdvZQoKYGBge3J9CmdvZV90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KGdvZV90ZXh0JHRleHQpLCAyKQoKcHJpbnQoZ29lX3RleHQpCmBgYAoKYGBge3J9CmVtb3Rpb25zX21ha2VyKGdvZV90ZXh0JHRleHQsICJnb2UiKQpgYGAKCiMjIyBoYXJ2YXJkCgpgYGB7cn0KaGFydmFyZF90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KGhhcnZhcmRfdGV4dCR0ZXh0KSwgMikKCnByaW50KGhhcnZhcmRfdGV4dCkKYGBgCgpgYGB7cn0KZW1vdGlvbnNfbWFrZXIoaGFydmFyZF90ZXh0JHRleHQsICJoYXJ2YXJkIikKYGBgCgojIyMgbGVpY2VzdGVyCgpgYGB7cn0KbGVpY2VzdGVyX3RleHQkc2VudGltZW50IDwtIHJvdW5kKGdldF9zZW50aW1lbnQobGVpY2VzdGVyX3RleHQkdGV4dCksIDIpCgpwcmludChsZWljZXN0ZXJfdGV4dCkKYGBgCgpgYGB7cn0KZW1vdGlvbnNfbWFrZXIobGVpY2VzdGVyX3RleHQkdGV4dCwgImxlaWNlc3RlciIpCmBgYAoKIyMjIG1hbmNoZXN0ZXIKCmBgYHtyfQptYW5jaGVzdGVyX3RleHQkc2VudGltZW50IDwtIHJvdW5kKGdldF9zZW50aW1lbnQobWFuY2hlc3Rlcl90ZXh0JHRleHQpLCAyKQoKcHJpbnQobWFuY2hlc3Rlcl90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcihtYW5jaGVzdGVyX3RleHQkdGV4dCwgIm1hbmNoZXN0ZXIiKQpgYGAKCiMjIyBtaXQKCmBgYHtyfQptaXRfdGV4dCRzZW50aW1lbnQgPC0gcm91bmQoZ2V0X3NlbnRpbWVudChtaXRfdGV4dCR0ZXh0KSwgMikKCnByaW50KG1pdF90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcihtaXRfdGV4dCR0ZXh0LCAibWl0IikKYGBgCgojIyMgc2IKCmBgYHtyfQpzYl90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KHNiX3RleHQkdGV4dCksIDIpCgpwcmludChzYl90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcihzYl90ZXh0JHRleHQsICJzYiIpCmBgYAoKIyMjIHN0YW5mb3JkCgpgYGB7cn0Kc3RhbmZvcmRfdGV4dCRzZW50aW1lbnQgPC0gcm91bmQoZ2V0X3NlbnRpbWVudChzdGFuZm9yZF90ZXh0JHRleHQpLCAyKQoKcHJpbnQoc3RhbmZvcmRfdGV4dCkKYGBgCgpgYGB7cn0KZW1vdGlvbnNfbWFrZXIoc3RhbmZvcmRfdGV4dCR0ZXh0LCAic3RhbmZvcmQiKQpgYGAKCiMjIyB0cmluaXR5CgpgYGB7cn0KdHJpbml0eV90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KHRyaW5pdHlfdGV4dCR0ZXh0KSwgMikKCnByaW50KHRyaW5pdHlfdGV4dCkKYGBgCgpgYGB7cn0KZW1vdGlvbnNfbWFrZXIodHJpbml0eV90ZXh0JHRleHQsICJ0cmluaXR5IikKYGBgCgojIyMgd3YKCmBgYHtyfQp3dl90ZXh0JHNlbnRpbWVudCA8LSByb3VuZChnZXRfc2VudGltZW50KHd2X3RleHQkdGV4dCksIDIpCgpwcmludCh3dl90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcih3dl90ZXh0JHRleHQsICJ3diIpCmBgYAoKIyMjIHlhbGUKCmBgYHtyfQp5YWxlX3RleHQkc2VudGltZW50IDwtIHJvdW5kKGdldF9zZW50aW1lbnQoeWFsZV90ZXh0JHRleHQpLCAyKQoKcHJpbnQoeWFsZV90ZXh0KQpgYGAKCmBgYHtyfQplbW90aW9uc19tYWtlcih5YWxlX3RleHQkdGV4dCwgInlhbGUiKQpgYGAKCiMjIyBDcmVhdGlvbiBvZiBhdmVyYWdlIHNlbnRpbWVudCBmb3IgZWFjaCBIRUkKCmBgYHtyfQojIENvbWJpbmUgYWxsIHRoZSBzZW50aW1lbnQgc2NvcmVzIGludG8gb25lIGRhdGEgZnJhbWUKYWxsX3NlbnRpbWVudHMgPC0gcmJpbmQoCiAgICBkdWtlX3RleHRbLCBjKCJpZCIsICJzZW50aW1lbnQiKV0sCiAgICBlcGZsX3RleHRbLCBjKCJpZCIsICJzZW50aW1lbnQiKV0sCiAgICBnb2VfdGV4dFssIGMoImlkIiwgInNlbnRpbWVudCIpXSwKICAgIGhhcnZhcmRfdGV4dFssIGMoImlkIiwgInNlbnRpbWVudCIpXSwKICAgIGxlaWNlc3Rlcl90ZXh0WywgYygiaWQiLCAic2VudGltZW50IildLAogICAgbWFuY2hlc3Rlcl90ZXh0WywgYygiaWQiLCAic2VudGltZW50IildLAogICAgbWl0X3RleHRbLCBjKCJpZCIsICJzZW50aW1lbnQiKV0sCiAgICBzYl90ZXh0WywgYygiaWQiLCAic2VudGltZW50IildLAogICAgc3RhbmZvcmRfdGV4dFssIGMoImlkIiwgInNlbnRpbWVudCIpXSwKICAgIHRyaW5pdHlfdGV4dFssIGMoImlkIiwgInNlbnRpbWVudCIpXSwKICAgIHd2X3RleHRbLCBjKCJpZCIsICJzZW50aW1lbnQiKV0sCiAgICB5YWxlX3RleHRbLCBjKCJpZCIsICJzZW50aW1lbnQiKV0KKQoKIyBDYWxjdWxhdGUgdGhlIGF2ZXJhZ2Ugc2VudGltZW50IGZvciBlYWNoIEhFSQpoZWlfYXZlcmFnZV9zZW50aW1lbnRzIDwtIGFnZ3JlZ2F0ZShzZW50aW1lbnQgfiBpZCwgZGF0YSA9IGFsbF9zZW50aW1lbnRzLCBGVU4gPSBtZWFuKSAlPiUKICAgIHJlbmFtZShhdmVyYWdlX3NlbnRpbWVudCA9IHNlbnRpbWVudCkKCiMgUHJpbnQgdGhlIHRhYmxlCnByaW50KGhlaV9hdmVyYWdlX3NlbnRpbWVudHMpCmBgYAoKIyBDbHVzdGVycwoKIyMjIEZ1bmN0aW9uIGZvciBjb3NpbmUgbWF0cml4CgpgYGB7cn0KY29zaW5lX21hdHJpeF9tYWtlciA8LSBmdW5jdGlvbih0YWJsZSl7CiAgbnVtZXJpY2FsX2NvbHVtbnMgPC0gdGFibGUgJT4lCiAgICBzZWxlY3QoLWlkKQogIAogICMgTm9ybWFsaXppbmcgY29sdW1ucwogIG5vcm1hbGl6ZWRfY29sdW1ucyA8LSBhcy5kYXRhLmZyYW1lKHNjYWxlKG51bWVyaWNhbF9jb2x1bW5zKSkKICAKICAjIFRyYW5zcG9zaW5nIHRoZSBkYXRhIHRvIGNvbXB1dGUgc2ltaWxhcml0eSBiZXR3ZWVuIGNvbHVtbnMKICB0cmFuc3Bvc2VkX2RhdGEgPC0gdChub3JtYWxpemVkX2NvbHVtbnMpCiAgCiAgIyBDb21wdXRpbmcgdGhlIGNvc2luZSBzaW1pbGFyaXR5IG1hdHJpeAogIHNpbWlsYXJpdHlfbWF0cml4IDwtIGFzLm1hdHJpeChwcm94eTo6ZGlzdCh0cmFuc3Bvc2VkX2RhdGEsIG1ldGhvZCA9ICJjb3NpbmUiKSkKICAKICAjIENvbnZlcnRpbmcgZGlzdGFuY2UgdG8gc2ltaWxhcml0eQogIHNpbWlsYXJpdHlfbWF0cml4IDwtIDEgLSBzaW1pbGFyaXR5X21hdHJpeAogIAogICMgUHJpbnQgdGhlIHNpbWlsYXJpdHkgbWF0cml4CiAgcHJpbnQoc2ltaWxhcml0eV9tYXRyaXgpCn0gIApgYGAKCiMjIyBDcmVhdGluZyB0YWJsZSBmb3IgY2x1c3RlciBhbGdvcml0aG1zCgpgYGB7cn0KIyBKb2luaW5nIGF0dHJpYnV0ZSBwZXJjZW50YWdlX3R3ZWV0cyAocGVyY2VudGFnZSBvZiB0d2VldHMgb3V0IG9mIGFsbCBwb3N0cykgYW5kIHBlcmNlbnRhZ2VfcmVwbGllcyAocGVyY2VudGFnZSBvZiByZXBsaWVzIG91dCBvZiBhbGwgcG9zdHMpIGZyb20gbnVtYmVyX3Bvc3RzIGFsc28gYWRkaW5nIHVuaXF1ZV9oYXNodGFncyAobnVtYmVyIG9mIHVuaXF1ZSBoYXNodGFncykgYW5kIGhhc2h0YWdfcGVyY2VudGFnZSAocGVyY2VudGFnZSBvZiBwb3N0cyB0aGF0IGNvbnRhaW4gYSBoYXNodGFnKSBmcm9tIGhhc2h0YWdzLCBwZXIgSEVJCmNsdXN0ZXJfdGFibGUgPC0gbWVyZ2Uoc2VsZWN0KGhhc2h0YWdzLCBpZCwgdW5pcXVlX2hhc2h0YWdzLCBoYXNodGFnX3BlcmNlbnRhZ2UpLCBzZWxlY3QoZGF0YV9yYXRpbywgaWQsIHBlcmNlbnRhZ2VfdHdlZXRzLCBwZXJjZW50YWdlX3JlcGxpZXMpLCBieSA9ICJpZCIsIGFsbD1UUlVFKQoKIyBKb2luaW5nIGF0dHJpYnV0ZSBhdmdfcG9zdHNfcGVyX2RheXMgKGF2ZXJhZ2Ugb2YgcG9zdHMgcGVyIGRheSkgZnJvbSBwb3N0c19wZXJfZGF5IHBlciBIRUkKY2x1c3Rlcl90YWJsZSA8LSBtZXJnZShjbHVzdGVyX3RhYmxlLCBzZWxlY3QocG9zdHNfcGVyX2RheSwgaWQsIGF2Z19wb3N0c19wZXJfZGF5cyksIGJ5ID0gImlkIiwgYWxsPVRSVUUpCgojIEpvaW5pbmcgYXR0cmlidXRlIGF2Z19wb3N0c19wZXJfd2Vla3MgKGF2ZXJhZ2Ugb2YgcG9zdHMgcGVyIHdlZWspIGZyb20gcG9zdHNfcGVyX3dlZWsgcGVyIEhFSQpjbHVzdGVyX3RhYmxlIDwtIG1lcmdlKGNsdXN0ZXJfdGFibGUsIHNlbGVjdChwb3N0c19wZXJfd2VlaywgaWQsIGF2Z19wb3N0c19wZXJfd2Vla3MpLCBieSA9ICJpZCIsIGFsbD1UUlVFKQoKIyBKb2luaW5nIGF0dHJpYnV0ZSBhdmdfcG9zdHNfaW5fYWNhZGVtaWNfdGltZSAoYXZlcmFnZSBvZiBwb3N0cyBkdXJpbmcgYWNhZGVtaWMgdGltZSkgIGZyb20gZGF0YV9wb3N0c19hY2FkZW1pYyBwZXIgSEVJCmNsdXN0ZXJfdGFibGUgPC0gbWVyZ2UoY2x1c3Rlcl90YWJsZSwgc2VsZWN0KGRhdGFfcG9zdHNfYWNhZGVtaWMsIGlkLCBhdmdfcG9zdHNfaW5fYWNhZGVtaWNfdGltZSksIGJ5ID0gImlkIiwgYWxsPVRSVUUpCgojIEpvaW5pbmcgYXR0cmlidXRlIGF2Z19wb3N0c19pbl92YWNhdGlvbl90aW1lIChhdmVyYWdlIG9mIHBvc3RzIGR1cmluZyB2YWNhdGlvbiB0aW1lKSBmcm9tIGRhdGFfcG9zdHNfdmFjYXRpb25zIHBlciBIRUkKY2x1c3Rlcl90YWJsZSA8LSBtZXJnZShjbHVzdGVyX3RhYmxlLCBzZWxlY3QoZGF0YV9wb3N0c192YWNhdGlvbnMsIGlkLCBhdmdfcG9zdHNfaW5fdmFjYXRpb25fdGltZSksIGJ5ID0gImlkIiwgYWxsPVRSVUUpCgojIEpvaW5pbmcgYXR0cmlidXRlIHRpbWVfb2ZfZGF5X3ZhbHVlIChudW1lcmljYWwgdmFsdWUgcmVmZXJyaW5nIHRvIHRpbWUgb2YgZGF5IHdoZXJlIGV2ZXJ5IEhFSSBtYWRlIG1vcmUgcG9zdHMpIGZyb20gZmF2b3VyaXRlX2hvdXJfaGVpIHBlciBIRUkKY2x1c3Rlcl90YWJsZSA8LSBtZXJnZShjbHVzdGVyX3RhYmxlLCBzZWxlY3QoZmF2b3VyaXRlX2hvdXJfaGVpLCBpZCwgdGltZV9vZl9kYXlfdmFsdWUpLCBieSA9ICJpZCIsIGFsbD1UUlVFKQoKIyBKb2luaW5nIGF0dHJpYnV0ZSB1cmxfcGVyY2VudGFnZSAocGVyY2VudGFnZSBvZiBwb3N0cyB0aGF0IGNvbnRhaW4gYW4gdXJsKSBmcm9tIHVybF91c2FnZSBwZXIgSEVJCmNsdXN0ZXJfdGFibGUgPC0gbWVyZ2UoY2x1c3Rlcl90YWJsZSwgc2VsZWN0KHVybF91c2FnZSwgaWQsIHVybF9wZXJjZW50YWdlKSwgYnkgPSAiaWQiLCBhbGw9VFJVRSkKCiMgSm9pbmluZyBhdHRyaWJ1dGUgYXZlcmFnZV9udW1fd29yZHMgKGF2ZXJhZ2UgbnVtYmVyIG9mIHdvcmRzIGluIHRoZSBwb3N0cykgZnJvbSBkYXRhX3Bvc3RzX2NvbnRlbnRfbWV0cmljcyBwZXIgSEVJCmNsdXN0ZXJfdGFibGUgPC0gbWVyZ2UoY2x1c3Rlcl90YWJsZSwgc2VsZWN0KGRhdGFfcG9zdHNfY29udGVudF9tZXRyaWNzLCBpZCwgYXZlcmFnZV9udW1fd29yZHMpLCBieSA9ICJpZCIsIGFsbD1UUlVFKQoKIyBKb2luaW5nIGF0dHJpYnV0ZSBhdmVyYWdlX3NlbnRpbWVudCAoYXZlcmFnZSBzZW50aW1lbnQgb2YgcG9zdHMpIGZyb20gaGVpX2F2ZXJhZ2Vfc2VudGltZW50cyBwZXIgSEVJCmNsdXN0ZXJfdGFibGUgPC0gbWVyZ2UoY2x1c3Rlcl90YWJsZSwgc2VsZWN0KGhlaV9hdmVyYWdlX3NlbnRpbWVudHMsIGlkLCBhdmVyYWdlX3NlbnRpbWVudCksIGJ5ID0gImlkIiwgYWxsPVRSVUUpCgpwcmludChjbHVzdGVyX3RhYmxlKQpgYGAKCmBgYHtyfQpjb3NpbmVfbWF0cml4X21ha2VyKGNsdXN0ZXJfdGFibGUpCmBgYAoKIyMjIEJhc2VkIG9uIHRoaXMgY29zaW5lIG1hdHJpeCB3ZSBkZWNpZGVkIHRvIHJlbW92ZSBhdmVyYWdlX251bV93b3JkcyAoaGlnaCBzaW1pbGFyaXR5IHRvIHVuaXF1ZV9oYXNodGFncywgaGFzaHRhZ19wZXJjZW50YWdlIGFuZCB1cmxfcGVyY2VudGFnZSksIHVuaXF1ZV9oYXNodGFncyAoaGlnaCBzaW1pbGFyaXR5IHRvIGhhc2h0YWdfcGVyY2VudGFnZSBhbmQgdGhlIGRpc3RhbmNlcyBvZiBoYXNodGFnX3BlcmNlbnRhZ2UgdG8gdGhlIG90aGVyIGNvbHVtbnMgaXMgbG93ZXIpIGFuZCBhdmdfcG9zdHNfcGVyX3dlZWtzLCBhdmdfcG9zdHNfaW5fYWNhZGVtaWNfdGltZSwgYXZnX3Bvc3RzX2luX3ZhY2F0aW9uX3RpbWUgYXJlIGFsc28gYWxsIHJlbW92ZWQgKGR1ZSB0byB0aGVpciBzaW1pbGFyaXR5IHRvIGF2Z19wb3N0c19wZXJfZGF5cykKCmBgYHtyfQpjbHVzdGVyX3RhYmxlIDwtIGNsdXN0ZXJfdGFibGUgJT4lCiAgc2VsZWN0KC1hdmVyYWdlX251bV93b3JkcywgLXVuaXF1ZV9oYXNodGFncywgLWF2Z19wb3N0c19wZXJfd2Vla3MsIC1hdmdfcG9zdHNfaW5fYWNhZGVtaWNfdGltZSwgLWF2Z19wb3N0c19pbl92YWNhdGlvbl90aW1lKQoKY29zaW5lX21hdHJpeF9tYWtlcihjbHVzdGVyX3RhYmxlKQpgYGAKCiMjIyBGdW5jdGlvbiBmb3IgY2x1c3RlciBtZXRob2QKCmBgYHtyfQpjbHVzdGVyX21ha2VyIDwtIGZ1bmN0aW9uKG51bV9jbHVzdGVycywgdGFibGUpewogICMgRXhjbHVkaW5nIGlkIGNvbHVtbiBmb3IgY2x1c3RlcmluZwogIGNsdXN0ZXJfZGF0YSA8LSBzZWxlY3QodGFibGUsIC1pZCkKICAgIAogICMgU2NhbGluZyB0aGUgZGF0YSBmb3Iga21lYW5zIG1ldGhvZAogIHNjYWxlZF9kYXRhIDwtIHNjYWxlKGNsdXN0ZXJfZGF0YSkKICAKICBrbWVhbnNfbW9kZWwgPC0ga21lYW5zKHNjYWxlZF9kYXRhLCBjZW50ZXJzID0gbnVtX2NsdXN0ZXJzLCBuc3RhcnQgPSAxMCkKCiAgIyBFeHRyYWN0aW5nIGNsdXN0ZXIgYXNzaWdubWVudHMKICBjbHVzdGVyX2Fzc2lnbm1lbnRzIDwtIGttZWFuc19tb2RlbCRjbHVzdGVyCiAgCiAgIyBDcmVhdGluZyBuZXcgZGF0YSBmcmFtZSBjb21iaW5pbmcgb3JpZ2luYWwgZGF0YSB3aXRoIGNsdXN0ZXIgYXNzaWdubWVudHMKICBjbHVzdGVyZWRfZGF0YSA8LSBjYmluZChjbHVzdGVyX3RhYmxlJGlkLCBjbHVzdGVyX2RhdGEsIGNsdXN0ZXIgPSBjbHVzdGVyX2Fzc2lnbm1lbnRzKQogIAogIGNsdXN0ZXJlZF9kYXRhIDwtIGNsdXN0ZXJlZF9kYXRhWywgYygiY2x1c3Rlcl90YWJsZSRpZCIsICJjbHVzdGVyIildCiAgCiAgcHJpbnQoY2x1c3RlcmVkX2RhdGEpCn0KYGBgCgojIyMgRnVuY3Rpb24gdG8gZGlzY292ZXIgYmVzdCBudW1iZXIgb2YgY2x1c3RlcnMKCmBgYHtyfQplbGJvd19tYWtlciA8LSBmdW5jdGlvbih0YWJsZSl7CiAgY2x1c3Rlcl9kYXRhIDwtIHNlbGVjdCh0YWJsZSwgLWlkKQogIHNjYWxlZF9kYXRhIDwtIHNjYWxlKGNsdXN0ZXJfZGF0YSkKICAKICB3c3MgPC0gdmVjdG9yKCkKICByYW5nZSA8LSAxOjEwCiAgCiAgZm9yIChrIGluIHJhbmdlKSB7CiAgICBrbWVhbnNfbW9kZWwgPC0ga21lYW5zKHNjYWxlZF9kYXRhLCBjZW50ZXJzID0gaywgbnN0YXJ0ID0gMTApCiAgICB3c3Nba10gPC0ga21lYW5zX21vZGVsJHRvdC53aXRoaW5zcwogIH0KICAKICBlbGJvd19kZiA8LSBkYXRhLmZyYW1lKGsgPSByYW5nZSwgV1NTID0gd3NzKQogIGdncGxvdChlbGJvd19kZiwgYWVzKHggPSBrLCB5ID0gV1NTKSkgKwogICAgZ2VvbV9saW5lKCkgKwogICAgZ2VvbV9wb2ludCgpICsKICAgIGxhYnMoeCA9ICJOdW1iZXIgb2YgQ2x1c3RlcnMiLCB5ID0gIldpdGhpbi1DbHVzdGVyIFN1bSBvZiBTcXVhcmVzIChXQ1NTKSIsCiAgICAgICAgIHRpdGxlID0gIkVsYm93IE1ldGhvZCBmb3IgT3B0aW1hbCBrIikgKwogICAgdGhlbWVfbWluaW1hbCgpCn0KYGBgCgojIyMgUGxvdCBvZiBFbGJvdyBNZXRob2QgYW5kIHNlbGVjdGlvbiBvZiBiZXN0IG51bWJlciBvZiBjbHVzdGVyIHRvIHZpZXcgaG93IEhFSXMgYXJlIGdyb3VwZWQKCmBgYHtyfQplbGJvd19tYWtlcihjbHVzdGVyX3RhYmxlKQpjbHVzdGVyX21ha2VyKDQsIGNsdXN0ZXJfdGFibGUpCmBgYA==